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ABSTRACT 


A  cognitive  radar  (CR)  system  is  one  that  observes  and  learns  from  the  environment,  then 
uses  a  dynamic  closed-loop  feedback  mechanism  to  adapt  the  illumination  wavefonn  so 
as  to  provide  system  perfonnance  improvements  over  traditional  radar  systems.  A  CR 
system  that  performs  multiple  hypothesis  target  classification  and  exploits  the  spectral 
sparsity  of  correlated  narrowband  target  responses  to  achieve  significant  performance 
improvements  over  traditional  radars  that  use  wideband  illumination  pulses  was  recently 
developed.  This  CR  system,  which  was  designed  for  Gaussian  target  responses,  is 
extended  to  non-Gaussian  targets. 

In  this  thesis,  the  CR  system  is  generalized  to  deal  effectively  with  arbitrary  non- 
Gaussian  distributed  target  responses  via  two  key  contributions:  (1)  an  important 
statistical  expected  value  operation  that  is  usually  evaluated  in  closed  form  is  evaluated 
numerically  using  an  ensemble  averaging  operation,  and  (2)  a  powerful  new  statistical 
sampling  algorithm  and  a  kernel  density  estimator  are  applied  to  draw  complex  target 
samples  from  target  distributions  specified  by  both  a  desired  power  spectral  density  and 
an  arbitrary  desired  probability  density  function.  Simulations  using  non-Gaussian  targets 
demonstrate  very  effective  algorithm  perfonnance.  As  expected,  this  perfonnance  gain  is 
realized  at  the  expense  of  increased  computational  complexity. 
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EXECUTIVE  SUMMARY 


Cognitive  radar  (CR)  is  a  concept  proposed  in  2006  by  S.  Haykin  to  improve  the 
performance  of  radar  systems  in  resource-constrained  and  interference-limited 
environments  [1].  A  CR  system  is  one  that  has  the  ability  to  observe  and  learn  from  the 
environment,  then  use  a  dynamic  closed-loop  feedback  mechanism  to  adapt  the 
illumination  waveform  so  as  to  provide  system  performance  improvements  over 
traditional  radar  systems.  Traditional  radar  systems  use  wideband  pulses,  or  chirps,  as 
illumination  waveforms.  These  have  an  advantage  in  that  they  ensure  that  all  of  the 
modes  of  the  target  are  excited.  The  disadvantage  is  that  these  wideband  illumination 
waveforms  use  more  energy  than  necessary  when  extended  target  responses  are 
correlated,  as  they  sometimes  are  in  practice.  Some  extended  targets  have  narrowband 
frequency  spectra;  much  of  the  spectrum  is  sparse.  Exciting  these  low-energy  parts  of  the 
spectrum  is  wasteful.  This  relative  spectral  sparsity  offers  opportunity  for  CR  design. 

Recent  research  by  R.  Romero  et  al.  developed  a  CR  system  that  exploits  the 
spectral  sparsity  of  correlated  narrowband  target  responses  via  real-time  adaptive 
illumination  waveforms  that  are  matched  to  the  targets  [2]— [6] .  These  matched 
waveforms  strategically  place  spectral  energy  in  the  parts  of  the  spectrum  that  correspond 
to  the  highest  energy  parts  of  the  target  spectra.  This  approach  provides  target 
classification  performance  advantages  for  the  CR  system.  For  individual  targets,  the  CR 
algorithms  of  Romero  et  al.  create  illumination  waveforms  that  are  optimal  in  the  sense 
that  they  maximize  the  signal-to-noise  ratio  (SNR)  in  the  radar  receiver.  In  a  multiple 
target  hypothesis  classification  setting,  the  individually  optimal  illumination  waveforms 
are  linearly  combined  in  a  weighted  sum  to  fonn  a  single  designed  illumination 
waveform  referred  to  as  a  probability  weighted  energy  (PWE)  waveform,  where  the 
weights  are  the  prior  probabilities  of  the  corresponding  classes.  Initially,  the  weights 
(priors)  are  assumed  equal.  An  iterative  maximum  a  posteriori  probability  (MAP) 
algorithm  is  used  to  update  the  priors  and  the  PWE  in  each  iteration.  In  practice,  after 
about  5  to  10  iterations,  the  updated  priors  converge  to  values  that  cause  the  posterior 


probability  density  of  the  measurements  given  the  target  to  be  maximized.  This  leads  to  a 
high  probability  of  correct  classification  (Pcc)  for  the  multiple  hypothesis  classification 
problem. 

The  CR  system  of  Romero  et  al.  is  very  effective.  It  was  derived  for  the  case  in 
which  both  the  target  responses  and  the  measurement  noise  are  Gaussian-distributed 
stochastic  processes.  For  the  purposes  of  this  thesis,  we  dub  this  system  the  Gaussian 
Classification  Cognitive  Radar  (GCCR)  system.  The  Gaussian  assumption  is  very 
appropriate  for  the  measurement  noise.  However,  in  some  cases,  it  may  not  be 
appropriate  for  the  target  response  signals.  Consequentially,  if  the  target  distributions 
vary  from  Gaussian,  the  GCCR  system  performance  may  be  compromised.  This  is 
especially  true  if  the  target  probability  densities  are  bi-modal  or  multi-modal.  There  exist 
strong  theoretical  and  operational  needs  in  the  CR  setting  for  the  ability  to  design 
adaptive  illumination  waveforms  for  non-Gaussian  targets. 

The  purpose  of  this  thesis  is  to  extend  the  CR  multiple  hypothesis  classification 
algorithm  of  Romero  et  al.  to  deal  effectively  with  arbitrary  non-Gaussian  distributed 
target  responses.  For  Gaussian  targets,  the  Bayesian  mathematics  of  the  GCCR  MAP 
algorithm  are  simplified  as  closed-fonn  expressions  for  the  posterior  density,  and  other 
terms  can  be  derived  with  some  effort.  For  arbitrary  non-Gaussian  targets,  the  basic 
Bayesian  mathematics  must  be  extended  appropriately  for  the  non-Gaussian  case,  then 
implemented  numerically  using  concepts  from  statistical  sampling  theory.  This  is  the 
primary  contribution  of  this  thesis.  One  key  result  of  this  research  is  the  recognition  that 
the  Bayesian  mathematics  can  be  extended  to  the  non-Gaussian  case  if  the  integrals  that 
evaluate  the  expected  value  operation  are  replaced  by  online  numerical  ensemble  average 
operations.  The  new  algorithm  derived  in  this  thesis  is  dubbed  the  Non-Gaussian 
Classification  Cognitive  Radar  (NGCCR)  system. 

Another  key  contribution  of  this  thesis  is  the  application  of  a  kernel  density 
estimator  and  a  new  statistical  sampling  algorithm  by  Nichols  et  al.  [7]  that  draws 
random  samples  from  a  desired  stochastic  process  that  is  specified  by  both  a  desired 
power  spectral  density  (PSD)  and  a  desired  arbitrary  (non-Gaussian)  probability  density 

function  (PDF).  The  sampling  algorithm  of  Nichols  et  al.  is  designed  for  real  data,  but  in 
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this  thesis  is  extended  for  complex  signals  by  applying  the  theory  of  complex  random 
variables.  Fundamental  statistical  sampling  algorithms  such  as  the  Metropolis-Hastings 
algorithm  and  the  Importance  Sampling  algorithm  draw  independent  identically- 
distributed  (IID)  samples  from  a  distribution  with  an  arbitrary  PDF.  However,  for  the  CR 
problem,  samples  from  a  distribution  that  is  specified  by  an  arbitrary  PDF  and  also 
correlated  by  the  PSD  of  the  narrowband  target  responses  are  required.  The  new 
algorithm  by  Nichols  et  al.  provides  the  capability  needed  for  this  work.  In  addition,  the 
sampling  algorithm  is  relatively  simple  to  implement,  making  it  very  valuable. 

The  new  CR  multiple  hypothesis  classification  algorithm  for  non-Gaussian 
distributed  target  classes  (the  NGCCR  system)  is  tested  using  simulated  data  sets  that  are 
non-Gaussian.  It  is  shown  that  the  Pcc  for  the  new  algorithm  using  the  MAP  designed 
PWE  illumination  waveform  is  significantly  greater  than  that  using  a  broadband  pulse 
waveform  over  a  wide  range  of  values  for  the  illumination  wavefonn  energy  Es .  As 

expected,  the  wideband  illumination  signal  and  the  PWE  give  approximately  the  same  Pcc 
when  the  transmit  energy  values  are  high  enough.  This  makes  intuitive  sense  and  is 
consistent  with  the  results  for  the  GCCR  system.  As  the  saying  goes,  “if  the  input  signal 
energy  is  large  enough,  then  any  wavefonn  is  a  good  waveform.”  In  practical  energy 
ranges,  the  new  NGCCR  system  provides  significant  perfonnance  improvement  over  the 
traditional  wideband  illumination  pulse,  making  it  valuable  in  a  CR  setting. 

The  tradeoff  for  the  new  NGCCR  system  is  that  it  involves  significantly  more 
computational  complexity  than  the  original  GCCR  system.  The  ability  to  deal  effectively 
with  non-Gaussian  target  responses  comes  at  the  price  of  more  computations.  This  is 
expected  and  necessary. 
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I.  INTRODUCTION  AND  MOTIVATION 


A.  BACKGROUND 

Traditional  radar  systems  use  simple  pulse  train  or  wideband  chirp  waveforms  as 
illumination  waveforms  to  obtain  target  responses.  The  task  of  classification  of  different 
targets  solely  relies  on  the  receiver  of  the  radar  system.  A  broadband  illumination 
waveform  can  cover  a  broad  frequency  spectrum  to  obtain  the  target  frequency  responses. 
However,  in  an  energy  constraint  environment,  if  the  target  set  is  auto-correlated  and  the 
spectra  are  narrowband,  much  illumination  energy  is  wasted  on  the  frequency  bands 
where  the  target  spectrum  is  small  or  zero.  In  the  cognitive  radar  setting,  the  radar 
transmits  signals  that  are  designed  to  match  the  characteristics  of  the  targets  for  improved 
target  recognition  results.  A  cognitive  radar  system  analyzes  the  statistics  of  the  dominant 
energies  of  an  ensemble  of  potential  targets  in  the  frequency  domain  and  formulates  the 
illumination  waveform  that  emphasizes  those  frequencies.  This  approach  can  improve  the 
efficiency  in  a  limited  energy  setting.  The  system  block  diagram  of  a  Cognitive  Radar 
(CR)  is  illustrated  in  Figure  1  [1]. 


Receiver 


Update  Prior 
Knowledge 


Statistical  Parameter  Estimates  and  Probabilistic 
Decisions  on  the  Targets 


Figure  1 .  Block  diagram  of  a  closed-loop  feedback  cognitive  radar  system. 

After  [1]. 
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Based  on  Figure  1,  a  cognitive  radar  system  consists  of  three  critical  ingredients: 
(1)  intelligent  signal  processing,  which  builds  on  learning  through  the  interaction  between 
the  radar  and  the  potential  targets;  (2)  feedback  from  the  receiver  to  the  transmitter, 
which  updates  the  intelligence;  and  (3)  preservation  of  the  information  content  of  the 
radar  returns  [1].  Cognitive  radar  distinguishes  itself  from  adaptive  radar  by  its  ability  to 
learn  and  adjust  to  the  environment.  Feedback  must  be  used  carefully  in  order  to  ensure 
stable  and  reliable  operation  at  all  times. 

The  design  of  illumination  waveforms  is  critical  to  the  performance  of  a  CR 
system.  The  illumination  signal  must  contain  enough  bandwidth  to  excite  the  resonant 
modes  of  the  target  in  order  to  obtain  a  full  characterization  of  the  target.  Moreover,  it 
must  contain  enough  power  to  obtain  useful  target  return  signals. 

B.  OBJECTIVE 

The  basic  approach  presented  by  R.  Romero  et  al.  [2]-[6]  designs  transmit 
adaptive  radar  illumination  waveforms  and  updates  the  prior  probabilities  of  the  targets 
for  classification  using  the  assumption  that  the  targets  are  Gaussian  distributed.  It  also 
solves  a  suite  of  automatic  target  recognition  problems,  using  both  the  Signal-to-Noise 
Ratio  (SNR)  Maximization  Method  and  the  Information-Theoretic  Method  [5].  This  work 
represents  the  current  state-of-the-art  in  the  area.  Here  we  extend  the  CR  system  to 
classify  non-Gaussian  distributed  targets.  This  research  is  focused  on  finding  a  general 
classification  algorithm  that  is  suitable  for  targets  with  any  statistical  distributions.  This 
research  requires  non-Gaussian  target  simulation,  probability  density  function  (PDF) 
estimation,  optimal  illumination  waveform  design  algorithms  based  on  SNR 
maximization,  and  target  classification  algorithms,  all  coded  using  MATLAB.  The 
research  aims  to  improve  the  probability  of  correct  classification  (Pcc)  with  transmit 
energy  constraints  for  the  non-Gaussian  target  scenario. 

C.  RELATED  WORK 

In  the  literature,  optimal  illumination  waveform  design  methods  have  been 

investigated  for  radar  systems  [2]— [6] .  Recent  research  by  R.  Romero  et  al.  [2]— [5] 

created  new  CR  algorithms  that  exploit  the  spectral  sparsity  of  correlated  narrowband 
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target  responses  to  design  illumination  waveforms  that  are  matched  to  the  targets.  These 
designed  waveforms  strategically  place  spectral  energy  in  the  parts  of  the  spectrum  that 
correspond  to  the  highest  energy  parts  of  the  target  spectra.  This  approach  provides  target 
classification  performance  advantages  for  the  CR  system.  For  individual  targets,  matched 
illumination  waveforms  are  optimal  in  the  sense  that  they  maximize  the  signal-to-noise 
ratio  (SNR)  in  the  radar  receiver  [5],  In  a  multiple  target  hypothesis  classification  setting, 
the  individually  optimal  illumination  waveforms  are  linearly  combined  in  a  weighted  sum 
to  form  a  single  designed  illumination  waveform  called  probability  weighted  energy 
(PWE),  where  the  weights  are  the  prior  probabilities  of  the  corresponding  classes. 
Initially,  the  weights  (priors)  are  assumed  equal.  An  iterative  maximum  a  posteriori 
probability  (MAP)  algorithm  is  used  to  update  the  priors  and  the  PWE  in  each  iteration. 
In  practice,  after  about  5  to  10  iterations,  the  updated  priors  converge  to  values  that  cause 
the  posterior  probability  density  of  the  measurements  given  the  target  to  be  maximized. 
This  leads  to  large  Pcc  for  the  multiple  hypothesis  classification  problem. 

The  CR  system  of  Romero  et  al.  [2]-[5]  is  very  effective.  It  was  derived  for  the 
case  in  which  both  the  target  responses  and  the  measurement  noise  are  Gaussian- 
distributed  stochastic  processes.  For  the  purposes  of  this  thesis,  we  dub  this  system  the 
Gaussian  Classification  Cognitive  Radar  (GCCR)  system.  The  Gaussian  assumption  is 
very  appropriate  for  the  measurement  noise.  However,  at  times,  it  may  not  be  appropriate 
for  the  target  response  signals.  Consequentially,  if  the  target  distributions  vary  from 
Gaussian,  the  GCCR  system  performance  may  be  compromised.  This  is  especially  true  if 
the  target  probability  densities  are  bi-modal  or  multi-modal.  There  exist  strong  theoretical 
and  operational  needs  in  the  CR  setting  for  the  ability  to  design  adaptive  illumination 
waveforms  used  for  non-Gaussian  targets. 

This  thesis  work  requires  the  ability  to  simulate  complex  target  signals  given 
desired  target  power  spectral  density  (PSD)  functions  and  PDF.  Previous  work  in  this 
area  lies  in  algorithms  for  Gaussian  stochastic  processes  and  algorithms  for  arbitrarily- 
distributed  (non-Gaussian)  stochastic  processes.  S.  Kay  [7]  presents  a  method  to  generate 
Gaussian  distributed  complex  signals  with  a  given  PSD.  This  technique  is  used  by 
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Romero  et  al.  in  the  GCCR  system  simulations.  For  the  work  in  this  thesis,  it  is 
necessary,  however,  to  generate  target  signals  with  desired  PSDs  and  desired  non- 
Gaussian  PDFs. 

A  key  contribution  of  this  thesis  is  the  application  of  a  new  statistical  sampling 
algorithm  by  Nichols  et  al.  [8]  that  draws  random  samples  from  a  desired  stochastic 
process  that  is  specified  by  both  a  desired  PSD  and  a  desired  non-Gaussian  PDF.  The 
sampling  algorithm  in  [8]  is  designed  for  real  data,  but  in  this  thesis  it  is  extended  for 
complex  signals  by  applying  the  theory  of  complex  random  variables.  Fundamental 
statistical  sampling  algorithms  such  as  the  Metropolis-Hastings  algorithm  and  the 
Importance  Sampling  algorithm  draw  independent,  identically  distributed  (IID)  samples 
from  a  distribution  with  an  arbitrary  PDF  [9],  [10].  However,  for  the  CR  problem,  it  is 
required  to  draw  samples  from  a  distribution  that  is  specified  by  an  arbitrary  PDF  and 
also  from  a  correlation  that  is  defined  by  the  PSD  of  the  narrowband  target  responses. 
The  new  algorithm  by  Nichols  et  al.  provides  the  capability  needed  for  this  work.  In 
addition,  the  sampling  algorithm  is  relatively  simple  to  implement,  making  it  very 
valuable. 

D.  APPROACH 

The  purpose  of  this  thesis  is  to  extend  the  CR  multiple  hypothesis  classification 
algorithm  of  Romero  et  al.  [2] — [5]  to  deal  effectively  with  arbitrary  non-Gaussian 
distributed  target  responses.  The  fundamental  classification  scheme  in  [3]  and  [4]  uses  a 
Bayesian  MAP  scheme  that  is  very  general,  in  the  sense  that  its  basic  formulation  is 
written  in  terms  of  PDFs  with  arbitrary  distributions.  However,  the  GCCR  approach  of 
[3]  and  [4]  specializes  the  Bayesian  formulation  by  assuming  Gaussian  target 
distributions.  This  specialization  allows  the  expected  value  to  be  written  as  an  integral  in 
closed  form.  This  integral  can  then  be  solved  in  closed  form  for  the  desired  posterior 
density. 

For  arbitrary  non-Gaussian  targets,  the  basic  Bayesian  equations  can  be  used  in 
their  most  general  form.  The  integrals  that  represent  the  expected  value  in  [3]  and  [4] 
cannot,  in  general,  be  solved  in  closed  form  as  they  were  in  the  GCCR  approach.  We 
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generalize  the  technique  by  implementing  the  expected  value  as  an  ensemble  average 
over  a  simulated  ensemble  of  target  realizations.  The  most  difficult  and  critical  problem 
lies  in  developing  an  algorithm  for  simulating  the  appropriate  target  samples  that 
constitute  the  ensemble.  This  problem  is  solved  and  described  in  the  thesis. 

The  key  contributions  of  this  research  he  in  using  the  most  general  Bayesian 
MAP  formulation,  in  which  the  integrals  that  evaluate  the  expected  value  are  replaced  by 
online  numerical  ensemble  average  operations.  The  new  algorithm  derived  in  this  thesis 
is  dubbed  the  Non-Gaussian  Classification  Cognitive  Radar  (NGCCR)  system. 

For  this  research,  all  simulations  are  coded  in  MATLAB,  and  all  MATLAB  m- 
files  and  functions  used  in  this  research  are  documented  in  the  Appendix. 

E.  THESIS  OVERVIEW  AND  ORGANIZATION 

The  matched  wavefonn  design  method  developed  by  R.  Romero  et  al.  for 
Gaussian  targets  (called  the  GCCR  system  in  this  thesis)  is  described  in  Chapter  II.  In  this 
chapter,  we  introduce  the  signal  and  system  model  for  the  radar,  define  the  power  spectral 
variance,  and  summarize  the  optimal  SNR  waveform  design  scheme  for  a  single  target. 

The  multiple  hypothesis  target  classification  method  used  in  the  GCCR  system  for 
Gaussian  targets  is  described  in  Chapter  III.  This  includes  closed  form  expressions  for 
posterior  densities  derived  using  the  Gaussian  assumption.  Then,  the  GCCR  classification 
theory  is  extended  to  deal  with  non-Gaussian  targets. 

The  motivation  and  mathematical  theory  of  statistical  target  signal  generation  is 
presented  in  Chapter  IV.  First,  the  concept  of  complex  signal  is  explained.  Next,  the 
algorithm  described  by  Kay  [7]  for  generating  a  complex  Gaussian  signal  with  a  specified 
PSD  is  presented.  Finally,  a  new  algorithm  by  Nichols  et  al.  is  presented.  It  provides  a 
relatively  simple  method  for  generating  a  real-valued  stochastic  process  with  a  specified 
PSD  and  a  specified  non-Gaussian  PDF.  We  then  show  that  this  algorithm  can  be 
extended  to  work  with  complex  stochastic  processes. 

The  specific  implementation  issues  involved  with  the  new  NGCCR  system  are 
presented  and  the  experimental  setup  used  to  test  the  new  classification  algorithm  is 
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described  in  Chapter  V.  Specifically,  we  describe  the  method  of  simulating  the  complex 
non-Gaussian  targets,  the  implementation  of  the  optimal  waveform  design  algorithm,  the 
MAP  multiple  hypothesis  classification  scheme,  and  the  algorithm  for  updating  the  prior 
probabilities  of  the  target  classes.  Finally,  the  Monte  Carlo  simulation  scheme  for  testing 
the  NGCCR  system  over  many  parameter  variations,  target  realizations  and  noise 
realizations  is  described. 

The  results  of  the  experiments  described  in  Chapter  V  for  the  NGCCR  system  are 
presented  in  Chapter  VI.  The  PSDs  and  PDF  estimates  of  the  simulated  targets  and 
representative  illumination  waveform  are  presented.  Finally,  curves  showing  the 
probability  of  correct  classification  vs.  illumination  signal  energy  for  both  the  optimal 
waveform  and  a  wideband  pulse  are  presented  and  compared.  This  curve  demonstrates 
the  perfonnance  improvement  offered  by  the  new  NGCCR  system. 

Conclusions  and  a  discussion  of  topics  for  future  research  are  presented  in 
Chapter  VII. 

The  Appendix  contains  the  MATLAB  m-files  and  functions  used  in  this  research. 
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II.  MATCHED  WAVEFORM  DESIGN 


In  this  chapter,  we  present  a  brief  summary  of  the  key  results  developed  by  Bell, 
Romero  et  al.  in  the  area  of  optimal  illumination  waveform  design  for  a  single  target  class 
[2]— [6] .  The  research  in  [2]-[6]  assumes  that  both  the  target  distribution  and  the  noise 
distribution  are  Gaussian.  For  this  thesis,  we  assume  that  the  targets  can  be  non-Gaussian 
or  arbitrarily  distributed,  but  the  noise  is  Gaussian. 

In  [5],  optimal  illumination  waveform  design  methods  were  developed  using  both 
mutual  infonnation  and  SNR  optimality  criteria.  Both  methods  are  effective.  In  this 
thesis,  we  focus  on  SNR  as  the  optimality  criterion  because  the  theoretical  derivations 
using  SNR  are  more  general  than  the  derivations  for  the  mutual  information  criterion. 
The  solutions  for  the  optimal  SNR  waveforms  in  [5]  are  valid  for  arbitrary  target 
distributions.  However,  the  derivations  using  mutual  information  require  the  assumption 
of  Gaussian  targets  [5].  Because  our  goal  is  to  use  arbitrarily-distributed  targets,  we  use 
the  SNR  approach. 

A.  SIGNAL  AND  SYSTEM  MODEL 

The  theoretical  radar  target  and  measurement  model  is  illustrated  in  Figure  2  [3]. 
The  complex  stochastic  radar  target  impulse  response  g(t )  is,  at  times  in  practice,  a 
finite-length  signal.  This  g(t )  can  be  modeled  as  a  stochastic  process  multiplied  by  a 
rectangular  window  of  duration  T  [3],  We  call  the  length  Ta  of  the  window  (in  seconds) 
the  mathematical  support  of  g(t) .  The  resulting  finite-length  impulse  response  is  called  a 
complex  stochastic  extended  target  signal  g(t) .  It  is  called  “extended,”  because  it  has 
finite  temporal  extent.  A  point  target  would  have  zero  extent  (an  impulse).  The  complex 
target  illumination  waveform  x(t)  with  duration  T  is  modeled  as  the  input  to  the  target 
impulse  response.  The  measured  complex  radar  return  signal  y(t )  is  modeled  as  the 
convolution  of  the  input  with  the  target  impulse  response  plus  additive  white  Gaussian 
noise  (AWGN)  n(t) 


y(f)  =  x(t)*g(t)  +  n(t). 
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Figure  2.  Complex-valued  baseband  signal  and  system  model  with  a  stochastic 

target.  After  [3]. 


We  note  that  in  [2]-[6],  the  concept  of  power  spectral  variance  (PSV)  rather  than 
PSD  is  used  to  describe  the  spectral  content  of  the  target  stochastic  processes.  This  is 
done  as  a  fonnal  way  to  deal  with  the  fact  that  we  have  finite-length  extended  target 
signals  that  cannot  strictly  be  considered  to  be  statistically  stationary.  When  the  PSD  is 
used  in  the  literature,  it  generally  is  defined  over  all  time  for  stationary  signals.  Because 
we  have  finite-duration  signals,  we  accommodate  them  using  the  concept  of  PSV.  The 
PSV  is  defined  in  terms  of  the  energy  spectral  density  (ESD)  given  by 


&(/>4|g,(/)| 


(2) 


Related  to  this  is  the  energy  spectral  variance  (ESV)  defined  as 


<(f)  =  E 


where  the  spectral  mean  juG(f)  of  the  target  process  is  defined  by 


(3) 


ft(/)  =  £[G, .(/)].  (4) 

The  key  point  of  this  approach  is  to  assume  that  the  finite-duration  target 
processes  are  stationary  within  their  mathematical  support  of  duration  To .  The  PSV  of  a 

target  signal  g(t)  is  then  defined  as  the  average  power  of  the  signal  over  the  time 
duration  Ta 
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(5) 


W)  = 


In  the  remainder  of  this  thesis,  the  PSV  of  a  signal  refers  to  this  definition. 


B.  OPTIMAL  SNR  WAVEFORM  DESIGN  FOR  A  SINGLE  TARGET 


Methods  for  designing  the  optimal  radar  illumination  waveform  for  a  single 
stochastic  target  using  the  SNR  optimality  criterion  are  described  in  [2]— [6] .  The  optimal 
SNR  waveform  x(t)  is  given  as  the  eigenvector  corresponding  to  the  maximum 
eigenvalue  of  the  eigenfunction  [3] 

4nax*(0  =  x(r)Rg  (( -  T)d T  (6) 

where  the  kernel  Rg (t)  is 

V')“f  crU/y^df.  (7) 

S  A[  J-  CO 

iV0 

The  PSD  of  the  AWGN  is  N0 ,  and  T  is  the  duration  of  the  illumination  waveform.  The 
optimal  illumination  waveform  for  a  single  target  is  denoted 

xop‘(t)  I  x(t).  (8) 

In  this  chapter,  we  reviewed  the  signal  and  system  model  for  radar  and 
summarized  the  design  method  for  the  optimal  SNR  illumination  waveform  for  a  single 
target  class.  In  the  next  chapter,  we  first  describe  how  this  optimal  waveform  for  a  single 
target  is  used  as  part  of  the  GCCR  multiple  hypothesis  target  classification  scheme, 
which  includes  a  prior  probability  update  that  is  developed  based  on  the  assumption  that 
the  targets  are  Gaussian  distributed.  We  then  introduce  a  new  generalized  target 
classification  scheme  that  does  not  rely  on  the  assumption  of  Gaussian  targets. 
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III.  MULTIPLE  HYPOTHESIS  TARGET  CLASSIFICATION 


In  this  chapter,  we  describe  the  MAP  multiple  hypothesis  [11  ]— [  13]  target 
classification  algorithms  under  consideration  in  this  thesis.  First,  we  describe  the  method 
used  in  the  GCCR  system,  which  was  designed  for  target  responses  that  are  Gaussian 
distributed  [2]— [5],  Then  we  describe  the  new  NGCCR  system,  which  relaxes  the 
Gaussian  assumption  and  is  shown  to  be  effective  for  arbitrary  target  distributions. 

We  assume  that  we  have  available  a  set  of  M  known  representative  target 
responses  |g;( t)}**  that  have  been  measured  or  simulated  in  advance.  Also,  assume  that 

we  calculate  in  advance  the  M  optimal  illumination  wavefonns  j associated 

with  the  M  target  classes.  We  then  use  radar  measurement  signals 
y(t)  -  x(t)*g(t)  +  n(t)  in  an  iterative  scheme  to  produce  an  illumination  waveform  that 
provides  high  classification  performance  measured  by  probability  of  correct  classification 
{Pcc)-  In  this  chapter,  we  describe  the  iterative  classification  algorithms  for  both  the 
Gaussian  and  non-Gaussian  target  assumptions. 

A.  DISCRETE-TIME  MEASUREMENT  MODEL 

The  measurement  model  described  earlier  is  written  in  terms  of  continuous  time. 
In  real-world  applications,  we  need  a  discrete  time  vector  formulation  which  we  describe 
next.  For  a  given  target  hypothesis  Hi ,  we  can  write  the  measurement  equation  as 

y  =  Xg,  +  n  (9) 

where  g,  is  Lx  1 ,  y  is  L  x  1  and  n  is  L  x  1 .  The  input  illumination  signal  x{t)  can  be 
written  as  an  L  x  1  input  vector  x .  This  input  vector  can  then  be  used  to  construct  the 
circulant  LyxL  circulant  convolution  matrix  Xk : 
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**(1)  0  .  0 

xk{2)  x,(l)  \  -  0 

*44)  *44-0  •••  *40  0 

0  *44)  *44-0  •••  *40 

:  o  x,(4)  •••  *42) 

:  :  0  : 

0  0  ...  0  xk(Lx) 

For  the  classification  problem,  we  use  multiple  illuminations  indexed  by  k  .  For  this  case, 
we  write  the  measurement  for  the  k,b  illumination  as 

lk=XkSi+!h-  (H) 

This  measurement  model  is  used  in  the  classification  algorithms,  as  we  next  show. 

B.  PWE  TRANSMISSION  TECHNIQUES 

In  a  multiple  target  hypothesis  classification  setting,  the  individually  optimal 
illumination  waveforms  {g,.(0}  are  linearly  combined  in  a  weighted  sum  to  form  a 

single  designed  illumination  waveform  called  PWE,  where  the  weights  are  the  prior 
probabilities  (or  priors)  Pt  of  the  M corresponding  hypotheses  Pt  =P(H ,)  [4],  For  each 

iteration,  a  new  illumination  wavefonn  is  transmitted  to  illuminate  the  targets.  The  priors 
and  the  PWE  are  updated  at  each  iteration  k  until  the  priors  converge  to  values  that 
optimize  the  posterior  density  that  are  described  later  in  this  chapter.  In  practice,  the 
algorithm  converges  to  useful  results  in  about  5  to  10  iterations.  The  desired  number  of 
iterations  is  chosen  in  advance  by  the  user,  so  no  additional  stopping  criterion  is  used.  In 
this  scheme,  the  transmitter  tends  to  place  the  illumination  waveform  energy  in  the 
spectral  bands  where  the  true  target  most  likely  resides.  This  multiple  transmission 
technique  is  tenned  the  PWE  technique  [4]  and  is  illustrated  in  Figure  3. 
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Figure  3.  PWE  multiple  transmission  technique  for  one  iteration.  From  [4]. 


The  PWE  multiple  transmission  technique  forms  individual  optimal  wavefonn 
spectra  |Tf(o;"(/)|  ,  which  are  developed  from  the  SNR  maximization  for  each  target 
hypothesis.  Then  each  individual  waveform  is  weighted  by  the  prior  probability  updates 
Rk  in  order  to  meet  an  overall  energy  constraint  Es  .  The  PWE  magnitude  spectrum 
squared  is  given  by 


M 


PWEi(f)  =  E,YJP,1  X°P'(f) 


2=1 


(12) 


where  P.k  is  the  prior  probability  of  the  ith  target  hypothesis  for  the  kth  transmission. 

This  PWE  illumination  waveform  is  the  sum  of  all  individual  matched  waveforms  scaled 
by  their  corresponding  updated  prior  probabilities.  In  the  case  where  the  individual 
matched  waveforms  are  designed  in  the  time  domain  using  an  eigenfunction  solution  (see 
Eq.  (6))  [3],  the  time  domain  PWE  waveform  is 


M  - 

pwEi(t)  =  jE;'£4kfxr(t). 


(13) 


2  =  1 

Note  that  the  PWE  is  scaled  by  a  square  root.  This  is  done  to  ensure  the  scaling  factor  for 
the  frequency  domain  has  the  same  scaling  effect  in  the  time  domain. 
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Note  that  the  PWE  is  used  as  the  radar  illumination  signal  to  illuminate  all  of  the 
targets.  Mathematically,  this  is  accomplished  in  the  measurement  model  by  writing  the 
PWE  signal  in  Eq.  (13)  as  an  L  x  1  input  PWE  vector  x  .  This  input  PWE  vector  is  used 
to  construct  the  circulant  convolution  matrix  Xk  as  described  in  Section  A. 

C.  ITERATIVE  ALGORITHMS  FOR  UPDATING  THE  PRIOR 

PROBABILITIES 

After  defining  some  key  quantities  and  relationships  regarding  the  measurements, 
we  now  review  the  MAP  algorithm  proposed  by  Romero  et  al.  [3]  for  updating  the  priors 
Pt  under  the  Gaussian  target  assumption  in  the  GCCR  system.  This  involves  using  closed 
form  expressions  for  the  required  PDFs  and  using  integrals  to  evaluate  the  required 
expected  value  operation.  Second,  we  derive  an  extension  of  this  MAP  algorithm  under 
the  assumption  that  the  target  random  variables  have  arbitrary  (non-Gaussian) 
distributions  for  the  NGCCR  system.  This  involves  replacing  the  integrals  in  the  GCCR 
system  with  a  numerical  ensemble  averaging  operation  to  evaluate  the  required  expected 
values.  It  also  requires  a  statistical  sampling  algorithm  to  draw  samples  from  the  desired 
target  distributions  specified  by  both  a  PDF  and  a  PSD.  This  sampling  algorithm  is 
described  in  a  succeeding  chapter. 

In  the  closed-loop  radar  multiple  hypothesis  classification  scheme,  for  the  ith 
target  for  (k  +  Y)th  transmission,  the  normalized  prior  probability  update  rule  is  given 
by  [3] 

pk+ 1  pki(yx,y2’-’yk)po 

i  M  (14) 

i= 1 

where  //(yl,y2,...,y/.)  is  the  joint  PDF  of  the  multiple  measurements  at  the  k'h 
illumination.  Initially,  the  weights  (priors)  are  assumed  equal  because  we  have  no  prior 
information  about  them.  We  then  let  the  iterative  algorithm  adjust  the  priors  appropriately 
according  to  the  MAP  criterion  [4],  [12],  [13].  The  initial  prior  P0  is  given  by 
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p»=17  (15) 

M 

where  M  is  the  number  of  the  target  hypotheses.  Next,  we  derive  the  prior  probability 
update  algorithms  required  in  the  CRGC  and  NGCCR  systems. 

1.  True  Measurements  vs.  Simulated  or  Estimated  Measurements 

Let  us  first  distinguish  clearly  between  true  measurements  and  simulated  (or 
estimated)  measurements .  By  a  true  measurement,  we  mean  the  measurement  that  a  radar 
system  would  acquire  in  an  operational  scenario.  We  denote  the  model  for  this  by 
yk  =  Xkgi  +  nk .  This  means  that  the  true  target  is  gt ,  and  we  received  a  measurement  of 

its  response  to  a  radar  illumination  signal. 

Let  us  now  consider  the  fact  that  for  our  problem,  we  have  prior  knowledge  of 

representative  target  responses  {  £,-(/)}  that  we  have  stored  in  memory.  We  can  later 

draw  samples  from  their  distributions  and  generate  an  ensemble  of  representative  signals 
that  we  can  use  in  our  classification  scheme.  We  can  simulate  “estimated  measurements” 
yk  =  Xkg:  that  are  used  as  an  estimate  of  the  measurement  one  would  expect  from  the 

known  target  classes.  Here,  we  let  the  “hat”  on  the  target  vector  gt  denote  the  concept  of 
a  single  simulated  realization  of  the  complex  random  target  vector  gt .  We  show  that  we 
can  use  these  estimated  measurements  in  a  MAP  scheme  to  update  the  target  priors. 

For  example,  consider  that  we  have  four  target  hypotheses  (M  =  4)  ,  and  the  third 
target  g:,  is  the  true  target  realization  corresponding  to  our  current  measurement 

yk  =  Xkg3  +  nk .  We  can  generate  four  estimated  measurements  |  v, =  Xkg^  |  that  can  be 
used  for  comparison  with  the  radar  measurement  yk  =  Xkg3  +  n k .  Next,  we  show  how  this 
comparison  is  done. 
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2. 


The  MAP  Scheme  for  Updating  the  Priors  in  the  GCCR  System 


The  iterative  update  algorithm  for  the  prior  probabilities  requires  that  we 
illuminate  the  target  multiple  times  and  measure  the  resulting  set  of  target  responses 
,  k  -  1,2, K  .  The  Bayesian  MAP  approach  requires  prior  knowledge  of 

l  1  m 

the  prior  PDFs  ,  where  M  is  the  number  or  target  classes.  These  priors  are 

available.  We  also  assume  that  the  M  target  responses  are  known  in  advance  because  we 
have  measured  them  or  modeled  them.  For  the  GCCR  system,  the  PDFs  have  been 
modeled  as  zero  mean  Gaussian  with  target  covariance  matrix  Ka  ,  denoted 

gj  :  N[0,  Kg  J  .  For  the  NGCCR  system,  we  estimate  these  PDFs  as  described  in  a 
succeeding  chapter. 

We  also  require  the  joint  conditional  PDF  pf(yl,y2—,yk  I  g,)  for  the  k'h  iteration. 

This  is  known  as  the  posterior  density  of  the  k  measurements  {\\,y2,...,yk\  conditioned 

on  a  particular  realization  of  the  ith  target  g; .  The  proper  formation  of  this  posterior  PDF 

is  the  key  to  the  effectiveness  of  the  MAP  algorithm.  The  k  measurements  conditioned 
on  the  target  realization  g(  are  independent,  so  their  joint  conditional  density  can  be 

written  as  the  product  of  the  individual  conditional  densities  p*  iyk  |g( )  as  follows: 


yk 


Si)=YlPi{yk\si). 

k-\ 


(16) 


The  key  to  the  GCCR  system  is  the  formulation  of  the  individual  posterior 
conditional  densities  | pf  (w  g,  j|  as  functions  of  the  true  measurements  \yk) k  ]  and 


the  estimated  measurements  j y.  k 


X k  g;  |  | .  Because  the  measurements  are  functions  of 


zero  mean  Gaussian  target  responses  and  zero  mean  noise  sequences,  Romero  et  al. 
consider  the  posterior  conditional  densities  as  Gaussian  PDFs  and  model  them  as  such 
[3]: 
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(17) 


Pt  [ik  |i, )  =  iL  |  exP  ~{yk  -  xdi  )H  kn  (: ik  -  x4 )  • 

71  L  J 

In  this  expression,  the  estimated  measurement  j yik  =  X  k gt j  ^  takes  on  the  role  of 
the  mean  of  the  Gaussian  random  variable  yk ,  and  the  covariance  matrix  K  v  takes  on  the 
role  of  the  covariance  Ky.  We  can  denote  this  as  follows:  yk  gt  □  N  [  Xkgj ,  K N  J  .  Note 

that  Kn  -  ■  I  where  Pn  is  the  noise  power,  and  /  is  an  identity  matrix  of  size 

Ly  x  Lv .  This  formulation  is  intuitively  appealing,  because  the  difference  yk  -  Xkg;  can 

be  interpreted  as  an  error  term  between  the  true  measurement  and  the  estimated 
measurement.  For  convenience,  call  this  error  term  %  k  =  yk  -  Xkg  .  Because  the 

measurement  model  is  yk  =  Xkg.t  +  nk  ,  if  g;  happens  to  be  the  true  target  in  the 
operational  scenario,  then  the  error  is  approximately  the  noise,  which  is  Gaussian 
distributed  n  :  N  [0,  K N  ] .  The  posterior  PDF  described  by  Eq.  (17)  models  this  noise 
very  nicely.  If  gt  happens  to  be  one  of  the  estimated  targets  (not  the  true  target)  in  the 
operational  scenario,  then  we  expect  the  error  to  be  somewhat  larger  than  it  is  when  gj  is 

the  true  target.  However,  the  error  can  be  thought  of  as  approximately  the  noise  with 
some  perturbation.  This  is  intuitively  appealing,  because  it  means  that  the  differences 
among  the  various  targets  manifest  as  varying  values  of  the  posterior  densities.  In 
addition,  Romero  et  al.  demonstrated  that  this  formulation  made  the  GCCR  system  very 
effective  for  multiple  hypothesis  target  classification  [3]— [5] . 

Continuing  with  the  GCCR  system  development,  we  can  write  the  posterior 
conditional  PDF  for  the  set  of  k  independent  measurements  as  [3] 

Pki[yl,y2,...,yk\g)  =  Y\ p-  (yk \gt ) 
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Given  the  quantities  described  above,  the  joint  density  of  the  multiple 
measurements  can  be  written  as  the  expected  value  of  the  joint  density  of  the  multiple 
measurements  conditioned  on  the  target: 


For  the  Gaussian  target  assumption,  we  can  write  the  expected  value  operation  in 
integral  fonn: 


>vA- 


g, 


)p(gi)dgl  • 


(20) 


Since  the  target  signal  is  based  on  a  Gaussian  assumption,  the  prior  density  of  the  target 
p(gi )  is  assumed  to  be  a  zero-mean  complex  Gaussian  signal,  and  its  PDF  is  given  by 


p(gi)  =  — 

1 

-exp 

H  T/~  —  \ 

-g,  Kg ,g, 

-  kl 

(21) 


The  final  expression  for  the  joint  PDF  is  derived  in  closed  form  by  evaluating  the 
expected  value  in  integral  form  as  written  in  Eq.  (20)  and  can  be  simplified  to  [3] 


Pi(yi,-,yk)  =  Eg,  Pi(y^-,yk  I &) 


0 


-i 


X 


K 


exp 


LK  jy-  K 

71  Kn 


K 


exp 


K 


-XziX'z* 


k=\ 


(22) 


k= 1 


K 


ZWz*  e~'2X^‘z<- 


k= i 


where  Kg  is  the  target  covariance  matrix  of  target  hypothesis  vector,  fij,  denotes  the 
expected  value  over  the  variable  g ,  and  Q  is  defined  by 


K 


Q=Kt+'£x“K»x>- 


(23) 


k=\ 
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The  expression  in  Eq.  (22)  is  used  in  the  prior  probability  update  algorithm  of  Eq.  (14). 
Once  the  priors  are  updated,  they  are  used  to  update  the  PWE  in  Eq.  (13).  This  PWE  is 
then  used  by  the  radar  as  the  new  illumination  signal  for  the  next  iteration. 


3.  Classification  Algorithm  for  the  NGCCR  System 


In  this  section,  we  derive  an  extension  of  the  MAP  algorithm  under  the 
assumption  that  the  target  random  variables  have  arbitrary  (non-Gaussian)  distributions 
for  the  NGCCR  system.  This  involves  replacing  the  integrals  in  the  GCCR  classification 
algorithm  with  a  numerical  ensemble  averaging  operation  to  evaluate  the  required 
expected  values. 


The  fundamental  Bayesian  MAP  approach  for  the  NGCCR  system  is  the  same  as 
that  for  the  GCCR  system.  We  desire  to  evaluate  the  joint  density  of  the  multiple 
measurements,  which  can  be  written  as  the  expected  value  of  the  joint  density  of  the 
multiple  measurements  conditioned  on  the  target: 


(24) 


Recall  that  for  the  GCCR  system,  this  can  be  written  in  closed  form  using  the 
expected  value  operation  in  integral  form: 

Pt  ili’li’  ■■■’lk)  =  \  Pi  ili’b’ ■■■’ Ik  | gt ) P(g,)dgr  (25) 

gi 


However,  for  the  NGCCR  system,  we  assume  that  the  target  responses  have 
arbitrary  (non-Gaussian)  distributions,  so  the  desired  PDFs  cannot  be  written  in  closed 
form.  As  a  result,  we  cannot  use  the  integral  form  of  the  expected  value  Eg  {•}  . 

Nonetheless,  the  posterior  joint  conditional  density  of  the  multiple  measurements  given  a 
realization  of  a  target  is  similarly  given  by  [3] 
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This  expression  remains  general  and  can  be  used  in  the  NGCCR  system  because  the 
targets  appear  as  the  variable  on  which  the  measurements  are  conditioned  and  the  prior 

densities  of  the  targets  p(gi )  do  not  appear  explicitly.  The  j  p(g,  )}.  |  do,  however,  come 

into  play  in  the  required  statistical  sampling  of  their  distributions  as  part  of  a  general 
complex  random  number  generator  described  in  a  succeeding  chapter.  The  multiple 
measurements  are  conditioned  on  single  realizations  of  the  target  vectors  g, ,  which  are 

treated  as  detenninistic  vectors. 

The  key  to  the  effectiveness  of  the  NGCCR  system  lies  in  the  numerical 
evaluation  of  the  expected  value  operation  E  {•{  in  Eq.  (24);  i.e.,  using  the  ensemble 

average  of  an  ensemble  of  many  realizations  of  the  posterior  joint  conditional  joint 
density  defined  by  Eq.  (26).  The  approach  is  to  create  Ng  realizations  of  the  target  g. 

and  use  these  to  form  an  ensemble  j/;C  [y],y2,...,yk  |g( )j  *  of  realizations  of  the  joint 

conditional  posterior  density.  We  can  then  compute  a  numerical  estimate  of  the  expected 
value  Eg  {•{  by  computing  the  ensemble  average  of  the  values  of  the  posterior  densities 
such  that 


Once  we  have  these  joint  densities,  we  use  them  to  update  the  prior  probabilities.  We  use 
the  updated  priors  to  update  the  PWE  illumination  signal,  which  is  then  transmitted  by 
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the  radar  to  illuminate  the  targets.  The  user  defines  in  advance  the  desired  number  of 
iterations  to  use.  When  that  number  of  iterations  is  complete,  the  final  target 
classification  decision  is  made. 

Numerically,  of  course,  this  evaluation  of  the  expected  value  comes  with  tradeoffs. 
It  allows  us  to  deal  effectively  with  non-Gaussian  distributed  target  classes.  This 
procedure  comes  at  the  price  of  increased  computational  complexity  and  depends  on  the 
number  of  target  realizations  generated  for  use  in  the  target  ensemble  above.  The 
purpose  of  this  thesis  is  to  show  the  derivation  of  the  classification  scheme  for  the 
NGCCR  system.  Future  research  should  focus  reducing  the  computational  complexity 
and/or  speed  of  implementation  of  the  NGCCR  system. 

In  this  chapter,  we  discussed  the  Gaussian  based  target  classification  scheme  used 
in  previous  researches  and  introduced  a  classification  method  that  eliminates  the 
Gaussian  assumption.  This  method  expands  the  application  of  the  target  classification 
scheme  for  complex  targets  with  any  arbitrary  density  distributions.  In  the  next  chapter, 
the  specific  implementation  of  the  experiment  is  outlined.  This  includes  the  discussion  of 
the  target  generation,  optimal  waveform  design,  and  target  classification  processes  used 
in  this  research. 
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IV.  TARGET  SIMULATION 


In  this  chapter,  we  introduce  the  concept  of  creating  stochastic  complex  signals. 
First,  the  feasibility  of  generating  meaningful  complex  signals  is  discussed.  Next,  we 
explain  the  theory  of  generating  Gaussian  complex  signals  with  a  specified  PSD.  Finally, 
the  method  of  generating  non-Gaussian  complex  signals  with  a  specified  PSD  that  is  used 
in  this  research  is  explained. 

A.  GENERATION  OF  STOCHASTIC  COMPLEX  SIGNALS 

1.  Motivation  for  Using  a  Complex  Data  Model 

In  radar,  sonar,  and  many  communication  systems,  complex  signals  are  used 
extensively  to  represent  the  in-phase  and  quadrature  components  of  the  demodulated  data. 
The  complex  envelope  represents  the  concatenation  of  the  in-phase  and  quadrature 
components  into  the  form  of  real  and  imaginary  quantities.  Using  the  complex  data  model 
simplifies  analysis.  This  method  is  similar  to  the  simplification  of  using  a  Fourier  series 
of  complex  exponentials  to  represent  sines  and  cosines. 

2.  Concept  of  Complex  Density 

The  goal  for  the  complex  density  is  to  express  the  PDF,  CDF,  and  all  related 
statistical  measures  as  a  function  of  the  complex  variables  [7].  The  known  results  for  real 
random  variables  and  those  obtained  from  complex  random  variables  must  match.  Let  the 
variable  x  be  a  2N  x  1  vector  of  real  and  imaginary  components  as  two  individual 
random  variables  defined  as 

x  =  [xR,Xj  f  (28) 

and  let  the  variable  x  be  a  single  complex  random  variable  defined  as 

X  =  XR+jxI.  (29) 
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If  the  covariance  matrix  Ki  of  the  original  real  2/Vx  1  vector  is  constrained  such 

that 

K  =K  =  K  /  2  (30) 

x,  R 

and 

KIR  =  -KTR1  =  KI!  2,  (31) 

then  the  duality  of  the  density  functions  for  the  real  random  variables  and  the  complex 
random  variables  can  be  written  as 

fk(u,v)=f:x(u+jv).  (32) 

The  first  moment,  or  the  mean  vector,  can  be  expressed  as 

E{x}U  E{xR}  + jE{xI}  =  m~.  (33) 

The  second  moment,  or  the  correlation  matrix,  can  be  expressed  as 

e{xxh}  =  R-x  =  RXr  +  RX1  +  j(RIR  - Rri )•  (34) 

Based  on  this  concept,  the  complex  signals  that  are  used  in  this  research  can  be 
analyzed  separately  by  splitting  the  signal  vector  into  its  real  and  imaginary  components 

m. 

B.  GENERATION  OF  COMPLEX  SIGNALS  WITH  SPECIFIED  SPECTRA 
AND  PDF 

1.  Generation  of  a  Complex  Gaussian  Signal  with  a  Specified  PSD 

In  previous  research,  the  target  classification  scheme  for  the  cognitive  radar 
system  was  focused  on  detecting  complex  Gaussian  distributed  signals  with  certain 
power  spectral  densities.  A  traditional  method  used  to  generate  this  type  of  signal  was 
developed  by  Steven  Kay  [7].  This  algorithm  is  based  on  the  concept  that  for  a  Gaussian 
random  process,  the  PSD  is  specified  by  the  covariance  matrix  R.  Since  this  covariance 
matrix  is  a  symmetric  Toeplitz  matrix,  the  eigenvalues  Ai  and  the  eigenvectors  vt  can  be 
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easily  found.  Let  Sxx(f )  be  the  PSD  of  x[n] ;  then  as  N  — »oo,  the  eigenvalues  A.  and  the 
eigenvectors  vi  are 

U .  =  5  (/) 

i  XX  v«/  i  7 

1  „  (35) 

L  =  — []  exp(y2^/:)  exp(y4^/:)  ...  exp(y2;r(A-lX/;)f 

for  z  =  0,1,...,  A  —  1  and  f  ]  -  /N-  Then  covariance  matrix  R  is  defined  as 

R  =  ZA/V,-V,^.  (36) 

i= 0 

If  the  data  set  {x[0],x[l],...,x[A-l]|  is  randomly  generated  based  on  a  Gaussian 

distribution  where  x[n] :  N(0,  cr2 ) ,  then  the  zero-mean  complex  Gaussian  discrete  target 
signal  g[ri\  with  a  specified  covariance  matrix  R  (or  PSD)  can  be  created  by 

g[n]  =  Vr  -x[n]  (37) 

for  n  =  0,l,—,N-  \ . 

Using  this  method,  we  can  generate  a  complex  Gaussian  distributed  sequence 
with  a  specified  PSD  to  simulate  the  target  signals  for  the  existing  target  detection 
scheme.  As  useful  as  this  method  is  for  simulating  complex  Gaussian  targets,  it  cannot  be 
used  to  generate  the  non-Gaussian  complex  target  signals  that  are  needed  for  this 
research.  Therefore,  a  technique  for  generating  complex  signals  with  both  specified 
arbitrary  PDF  and  PSD  is  required. 

2.  Generate  a  Complex  Signal  with  a  Specified  PDF  and  PSD 

Sampling  methods  such  as  the  Importance  Sampling,  Metropolis-Hasting 
Sampling,  and  the  Slice  Sampling  create  uncorrelated,  independent  and  identically 
distributed  (IID)  samples.  For  the  purpose  of  this  research,  a  correlated  non-Gaussian 
distributed  target  is  desired  to  evaluate  the  classification  system.  The  traditional  “Kay” 
method  of  generating  Gaussian  stochastic  signals  with  specified  PSDs  discussed  in  the 
previous  section  cannot  be  used  to  generate  non-Gaussian  distributed  signals.  Much 
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research  has  been  done  to  solve  the  problem  of  how  to  generate  observations 
x(n),n  =  0...N  - 1  with  a  given  PDF  and  auto-covariance  (or  PSD).  Although  the  details 
vary  in  different  methods,  all  approaches  tend  to  follow  a  similar  prescription  [8],  [14]— 
[17].  One  algorithm  developed  by  Nichols  et  al.  for  generating  spectrally  colored,  non- 
Gaussian  complex  signals  is  chosen  to  use  for  the  classification  system  [8],  However,  it 
must  be  noted  here  that  the  PDF  p(x )  and  the  PSD  Sxx  (/)  cannot  be  specified 
independently  because  these  two  properties  are  linked  through  the  mean  x  and  the 
variance  <j2x  of  the  signal. 

A  block  diagram  that  generally  summarizes  the  procedure  of  generating  signals 
with  both  a  desired  PDF  and  a  PSD  is  illustrated  in  Figure  4  [8].  First,  a  spectrally  white, 
Gaussian  distributed  sequence  is  linearly  filtered  by  the  specified  auto-covariance 
(equivalently  PSD)  to  obtain  a  sequence  with  a  Gaussian  PDF  and  the  specified  PSD. 
This  step  makes  use  of  the  fact  that  if  a  Gaussian  sequence  is  filtered  by  a  linear  system, 
the  output  is  a  Gaussian  sequence  as  well.  The  data  are  then  subject  to  a  zero-memory, 
nonlinear  (ZMNL)  transformation  in  order  to  produce  a  signal  with  the  desired  non- 
Gaussian  PDF. 


Linear  Digital  Filter 


ZMNL  Transformation 


x(n)  ~  N(0. 1)  auto  coirelation  (PSD)  PDF  and  PSD 


Figure  4.  General  process  for  generating  signals  with  a  user  defined  PDF  and  PSD. 

From  [8]. 


The  goal  of  the  algorithm  is  to  produce  a  sequence  of  observations  x[ri\ 
n  =  0,...,N-\  with  PDF  p(x)  and  a  PSD  Sxx(  f  ) .  First,  a  desired  PSD  function  Sxx(fk ) 
is  proposed  with  N  discrete  frequencies  at  the  bandwidth  fk  =(k- TV/ 2) A  ,  where 
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k  =  0,...,N- 1 .  The  frequency  bin  A/  is  dictated  by  the  temporal  sampling  interval  At 
such  that  Af  =  1/ (NAf) .  In  general,  Ar  should  be  chosen  in  accordance  with  the  Nyquist 

criterion  for  the  maximum  resolvable  frequency.  For  the  purpose  of  this  research,  the 
frequency  domain  is  normalized  to  1  Hz  for  the  analysis  of  all  signal  spectra.  Thus,  the 
frequency  bin  is  A  f  =  1  IN,  and  the  temporal  sampling  interval  is  At  —  1 .  The  discrete 
Fourier  transform  and  the  inverse  Fourier  transform  are  defined  as 


N- 1 

X[k]  =  FT(x[n])  = 

n= 0 


(38) 


and 


x[n]  =  FT~\X[k])  =  —YJX[k]ej2ukn,N  . 

X  k=o 


At  this  point,  the  relationship  between  the  signal  x[«]and  its  PSD  Sxx  (/)  in  tenns 
of  mean  and  variance  is  defined  by  [8] 


and 


-  JEM 

NAt 


(40) 


(41) 

1  k= 0 

A  data  sequence  s[n]  is  generated  from  a  desired  distribution  p(x) .  It  is  natural 
to  have  the  proposed  PDF  p(x)  to  specify  the  mean  of  the  signal  and  to  have  the 
proposed  PSD  Sxx(f )  to  specify  the  variance  of  the  signal.  Then  the  variance  of  the 
signal  s[n]  must  be  adjusted  to  conform  to  the  variance  of  the  proposed  PSD  Sxx  (/)  .  The 

variance  of  the  signal  is  a]  =  ^  ^  fy[n]-  T)2 .  With  the  variance  of  the  PSD  obtained 

from  Eq.  (41),  the  signal  s[n\  can  be  adjusted  to  a  scaled  sequence  x0[n\  that  conforms 
to  the  variance  of  Sxx  (/)  using 
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The  final  output  signal  needs  to  have  the  mean  specified  by  the  original  proposed 
PDF.  At  this  point,  the  mean  of  the  sequence  x0[n]  is  removed  for  that  purpose.  In  the 
last  step  of  the  procedure,  the  mean  is  put  back  into  the  signal.  The  final  step  for  the 
initialization  process  is  to  rearrange  the  values  in  x0[/J ,  n  =  0,...,N—  1  from  smallest  to 

largest  and  to  obtain  the  indexing  In  for  the  sequence  x0[n]  so  that 
x0[I0]  <  x0[/J  < ...  <  XotAv-i]  • 

After  the  initialization  process,  the  sequence  x0[n]  is  generated  with  the  target 
PDF  p(x)  and  the  variance  of  the  target  PSD  Sxx  (/) .  The  next  step  is  to  generate  a  new 
signal  x[n\  that  captures  the  phases  of  the  time-series  x0[n]  and  the  Fourier  amplitudes 
Xt[k]  of  the  proposed  PSD.  The  Fourier  amplitudes  of  the  PSD  Sxx  (/)  are  given  by 

X,m  =  JsJ7X T.  k  =  0,...,N-l,  (43) 

and  the  phases  of  the  values  in  the  time-series  x0[n]  are  given  by 


(f)\k]  =  tan 


f  Im(X0[k])^ 

Ue(X0[k])J  ’ 


(44) 


A  new  signal  x[n\  is  created  with  the  phases  of  xn[/?]  and  the  Fourier  amplitudes  of  the 
proposed  PSD  from 


i  N- 1 

*[»]=— E< 

k= 0 


Xt[k]e 


i7.nkn.IN 


FT~\eimX  [k]) . 


(45) 


This  signal  has  the  same  PSD  as  the  proposed  Sxx  (/) .  However,  due  to  the  central  limit 

theorem  [11],  the  process  of  the  inverse  Fourier  transfonn  has  the  effect  of  whitening  the 
signal  towards  a  Gaussian  distribution. 

At  this  point  in  the  process,  the  signal  x[n]  is  approximately  Gaussian  distributed 
and  has  the  desired  narrowband  PSD.  The  signal  x0  [/?]  has  the  target  PDF  and  a  constant 


PSD  (indicating  that  it  is  uncorrelated  or  white).  The  next  step  is  to  apply  a  nonlinear 
transfonnation  that  has  the  effect  of  auto-correlating  the  signal  x0  [/?]  and  ensuring  that  it 
has  approximately  the  desired  target  PSD  without  altering  its  PDF.  The  nonlinear 
transformation  involves  sorting  (or  rank  reordering)  the  values  in  x0  [/?]  ,  and  this 

reordering  process  causes  the  reordered  sequence  to  be  correlated.  Specifically,  the 
smallest  value  in  x0[n]  gets  assigned  the  same  position  of  the  smallest  value  in  x\n\ ,  and 

the  second  smallest  value  in  x0[n]  gets  assigned  the  same  position  of  the  second  smallest 
value  in  x[n\  etc.  The  sorted  values  of  x[n\,  from  smallest  to  largest,  are  denoted  x[Jn\ , 
where  x[/0]  <  x[J,  ]  < ...  <  x[Jv_, ] .  Then  the  signal 

x[J„]  =  x0[In],n  =  0,...,N-\  (46) 

by  definition  has  the  desired  target  PDF  and  approximately  the  desired  PSD.  This 
reordering  procedure  attempts  to  simulate  the  effect  of  the  ZMNL  transformation,  which 
is  a  process  of  mapping  a  signal  with  samples  drawn  from  the  target  PDF  into  another 
signal  that  possesses  the  target  PSD. 

Nichols  et  al.  developed  this  method  with  the  intention  of  creating  real  signals 
with  both  the  target  PDF  and  PSD.  However,  the  simulated  radar  target  signals  that  are 
used  in  this  research  must  be  complex  signals.  Based  on  the  concept  discussed  in  the  first 
section  of  this  chapter,  a  non-Gaussian  and  colored  spectrum  complex  signal  is  created  by 
combining  a  real  signal  and  an  imaginary  signal  such  that  both  are  generated  individually 
from  the  same  non-Gaussian  PDF  and  the  same  proposed  PSD.  With  the  combination  of 
these  two  concepts,  a  non-Gaussian,  colored  spectral  complex  signal  is  generated  to 
simulate  the  target  signals  in  this  research. 

For  the  work  done  in  this  thesis,  the  desired  specified  PSVs  are  symmetric  in  the 
frequency  domain.  This  was  done  because  of  a  limitation  of  the  statistical  sampling 
algorithm  of  Nichols  et  al.  The  algorithm  produces  only  real  sequences,  but  for  our  work, 
we  require  complex  sequences.  As  described  above,  we  form  a  complex  sequence  from 
two  real  sequences  generated  using  the  Nichols  algorithm.  The  PSVs  of  these  sequences 
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are  always  symmetric.  However,  for  some  operational  scenarios,  we  need  to  deal  with 
asymmetric  PSVs. 

Some  preliminary  work  was  done  to  create  asymmetric  PSVs  by  scaling  the  PSV 
peaks.  This  was  effective,  but  it  introduced  another  problem.  Once  the  scaled  PSV  is 
created,  the  only  way  to  use  it  to  create  a  time-domain  wavefonn  is  by  applying  an 
inverse  Fourier  transfonn  (IFT)  to  the  PSV.  This  provides  a  time-domain  waveform,  but 
it  is  not  rigorous  because  one  cannot  theoretically  compute  the  IFT  of  a  PSV.  The  PSV 
contains  only  the  modulus  of  the  complex  frequency  domain  quantity  with  no  phase 
information.  This  means  that  the  IFT  is  not  rigorously  defined.  Future  work  should  focus 
on  creating  a  generalized  method  for  sampling  a  distribution  specified  by  both  a  desired 
asymmetric  PSV  and  a  desired  PDF.  This  will  likely  involve  using  a  phase  retrieval 
algorithm  for  evaluating  the  IFT  of  a  PSV.  A  very  effective  phase  retrieval  algorithm  is 
described  in  [18]. 

The  concept  of  generating  complex  signals  needed  to  simulate  the  potential  target 
signals  for  this  research  was  explained  in  this  chapter.  Two  specific  methods  developed 
by  Kay  and  Nichols  respectively  were  examined.  This  research  uses  the  latter  method  to 
generate  complex  non-Gaussian  and  colored  spectral  signals  for  the  classification  system. 
In  the  next  chapter,  the  adaptive  waveform  design  method  used  in  this  research  is 
thoroughly  explained. 
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V.  EXPERIMENTAL  SETUP 


Experiments  and  tests  conducted  in  this  research  require  the  generation  of  non- 
Gaussian  colored  spectral  complex  target  signals,  optimal  wavefonn  computation  based 
on  proposed  target  PSVs,  coding  of  the  new  joint  posterior  PDF  estimate  algorithm,  and 
last  but  not  least,  a  Monte  Carlo  simulation  system  for  test  evaluations.  All  simulations 
are  conducted  in  MATLAB  and  are  explained  in  this  chapter. 

A.  TARGET  SIMULATION 


The  first  step  for  the  experiment  is  to  use  the  solution  in  [8]  to  define  a  number  of 
target  hypotheses  that  are  characterized  by  their  PSVs  as  experimental  subjects.  For  this 
research,  the  number  of  target  hypotheses  in  each  classification  test  set  is  set  to  four.  All 
proposed  PSVs  for  the  experiment  are  given  by  the  mixed-Gaussian  probability  density 
function 


S  (/)  = 


signal 


u-my 

2a1 


+  e 


(/+/')- 

2a1 


(47) 


where  Psignal  defines  the  total  power  in  the  illumination  waveform,  //  defines  the  band¬ 
pass  frequencies  of  the  two-sided  spectral  density  on  a  normalized  frequency  scale,  and 
a  defines  both  the  bandwidth  and  the  amplitudes  of  the  energy  bands.  It  must  be  noted 
here  that  the  bandwidths  and  the  amplitudes  of  the  energy  bands  in  each  PSV  are  related 
because  the  overall  power  constraint  for  the  PSV  is  set  by  Psi  al .  Thus,  changing  the 

parameter  o  either  increases  the  amplitude  and  narrows  the  energy  bands  or  decreases 
the  amplitude  and  expands  the  energy  bands.  Throughout  the  experiment,  the  target-to- 
noise  ratio  (TNR)  is  set  to  10  dB.  The  noise  in  the  receiver  is  set  to  unity  variance  with  a 
power  of  1  unit.  Thus,  P.  al  is  set  to  10  units  to  ensure  a  TNR  of  10  dB.  The  parameters 


u  and  a  are  the  two  tuning  variables  for  creating  different  PSVs.  An  example  of  a 
proposed  PSV  is  illustrated  in  Figure  5.  In  this  particular  example,  the  PSV  is  defined  by 
N  =  1024  data  samples  and  the  two  parameters  ju  =  0.2  and  cr  =  0.03 . 
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Figure  5.  An  example  of  a  PSV  defined  by  a  mixed  Gaussian  distribution 

Using  a  PSV  such  as  the  one  illustrated  in  Figure  5  along  with  the  MATLAB 
random  number  generator  based  on  a  chosen  density  distribution,  we  used  the  algorithm 
created  by  Nichols  et  al.  described  in  Chapter  IV  to  generate  real  data  samples  from  the 
proposed  distribution  with  a  PSV  that  is  very  similar  to  the  proposed  PSV.  In  the  case  for 
generating  one  complex  target  signal,  the  target  generator  combines  one  set  of  data 
samples  generated  for  the  real  component  with  another  set  of  data  samples  generated  for 
the  imaginary  component.  Taking  into  consideration  computational  speed,  we  set  the 
number  of  samples  used  through  the  experiment  to  N  =  64 .  The  last  step  for  target 
generation  is  to  scale  the  data  samples  properly  to  ensure  that  the  signal  conforms  to  the 
signal  power  defined  by  Psi  , .  Since  complex  signals  split  the  energy  between  the  real 

part  and  the  imaginary  part  and  the  scale  factor  in  the  frequency  domain  is  the  square  root 
of  the  scale  factor  in  the  time  domain,  the  proper  scaling  for  complex  signals  is  given  by 
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(48) 


g scaled  ^  ^ 


g(0 

4in' 


B.  SNR  MATCHED  WAVEFORMS  DESIGN 

As  soon  as  the  four  different  target  PSVs  are  defined,  the  optimal  individual 
waveforms  for  each  of  the  four  particular  target  hypotheses  can  be  generated.  The 
matched  waveform  design  algorithm  was  explained  in  Chapter  II.  To  implement  this 
SNR  maximization  algorithm,  first  the  covariance  matrix  R  of  each  PSV  Sxx(ft )  must  be 
computed.  For  AWGN,  the  expression  is  given  by  Eq.  (7).  An  alternative  way  to  compute 
the  integral  is  by  taking  the  complex  eigen  decomposition 

R  =  (49) 

i=0 

where  the  eigenvector  v;  is  given  by 

K-  =  exP(J2^(0/,)r  (50) 

for  i  =  0 ,...,  A -1  and  fi=~^  ■  Here,  f]  is  nonnalized  from  0  Hz  to  1  Hz. 

The  matched  waveform  for  a  particular  target  hypothesis  based  on  its  PSV  Sxx{f ) 
is  just  the  eigenvector  that  corresponds  to  the  maximum  eigenvalue,  which  is  the 
maximum  value  in  the  PSV  Sxx(f ) .  Here,  it  is  important  to  clarify  that  the  resulting 
optimal  waveform  generated  from  this  design  algorithm  is  a  temporal  waveform  because 
in  the  next  step,  the  probability  update  rule  is  different  for  waveforms  that  are  in  different 
domains. 

Last  but  not  least,  the  total  illumination  wavefonn  is  created  by  summing  all 
individual  optimal  waveforms  that  are  derived  from  each  of  the  four  particular  target 
hypothesis.  The  illumination  waveform  is  updated  by  a  probability  update  rule  for  each 
transmission.  However,  for  the  initial  transmission,  the  illumination  wavefonn  weighs 
each  target  hypothesis  equally.  The  illumination  waveform  is  given  by 
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i=  1 


where  P0-  %  defined  by  Eq.  (15). 

C.  TARGET  CLASSIFICATION  SCHEME 

After  the  illumination  waveform  is  generated  (along  with  a  target  realization  g 

from  an  unknown  target  hypothesis),  the  next  step  is  to  simulate  the  measurement  signal 
at  the  radar  receiver  by 

Z measured  ^XSi+H,  (52) 

where  n(t)  is  the  complex  white  noise  in  the  receiver. 

The  goal  for  the  experiment  is  to  test  the  new  classification  scheme.  This 
classification  algorithm  is  an  estimate  of  the  MAP  detection  scheme,  the  details  of  which 
were  explained  in  Chapter  III.  To  implement  this  algorithm,  simulated  target  responses 
from  each  target  hypothesis  are  needed  and  are  given  by 

y,=x%!  (53) 

for  i  =  l,. ..,4  and  gi  is  a  target  realization  based  on  one  particular  target  hypothesis. 
Then  each  one  of  the  four  joint  PDFs  of  the  measurement  y  for  each  target  is  estimated 
by  Eq.  (27).  The  covariance  matrix  KN  of  the  noise  is  given  by 


where  /  is  an  identity  matrix  of  size  Lv  x  L  .  Because  the  MAP  detection  scheme  makes 

the  decision  based  on  the  relative  probabilities  computed  from  the  conditional  PDFs,  all 
constants  in  the  expression  that  are  shared  by  the  four  joint  PDF  computations  can  be 
disregarded.  Thus,  for  the  k'1'  transmission  and  ith  hypothesis,  the  expression  in  Eq.  (27) 
can  be  simplified  to 
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(55) 


1  J 

. ^)  =  7ZexP 

J  7=1 


a: 


^iyk-yifKNiyk-yj) 


k= 1 


where  J  is  a  user  defined  number  that  specifies  the  number  of  realizations  in  the  target 
ensemble  and  yj  is  the  simulated  target  response  based  on  the  j'h  realization  of  the  i"' 
target  hypothesis  gt(t).  This  is  numerical  implementation  of  the  expected  value  using 

ensemble  averaging.  It  must  be  noted  here  that  the  parameter  J  dictates  the  accuracy  of 
the  expected  value  estimate,  which  means  that  a  higher  value  of  J  yields  a  more  accurate 
estimate  of  pk.(yv...,y,  g. ) .  However,  increasing  J  also  increases  the  computational 
time  for  each  transmission,  which  is  undesirable  for  target  detection  in  a  real  world 
scenario.  For  this  experiment,  the  value  J  is  chosen  as  40  to  balance  the  tradeoffs 
between  accuracy  and  speed;  J  was  found  empirically. 

The  last  step  for  the  target  classification  scheme  is  to  update  the  prior  probability 
for  each  target  hypothesis.  The  update  rule  is  given  by  Eq.  (14).  With  the  new  prior 
probability  weights  computed  for  each  target  hypothesis,  the  matched  waveform  for  the 
next  transmission  is  also  updated  by  Eq.  (13).  At  this  point,  one  radar  transmission  loop  is 
completed.  The  pseudo-code  block  diagram  of  one  radar  classification  evolution  is 
illustrated  in  Figure  6.  This  includes  the  transmission  of  the  first  matched  waveform, 
receiver  data  simulation,  target  classification,  and  the  probability  update  for  the  next 
transmission.  The  radar  simulator  performs  a  user  specified  number  of  transmissions  and 

makes  a  final  decision  for  target  classification  based  on  the  last  computed  prior  j/^j  . 

The  target  classification  decisions  are  then  used  to  compare  with  the  true  hypothesis  in 
order  to  evaluate  the  performance  of  the  radar  system.  This  concludes  one  cycle,  termed 
“noise  realization,”  for  the  target  classification  experiment. 


D.  MONTE  CARLO  SIMULATION  SYSTEM  SETUP 


Because  this  research  deals  with  stochastic  targets,  the  performance  of  the 
classification  system  cannot  be  truly  evaluated  from  one  target  realization.  A  Monte 

Carlo  simulation  is  set  up  for  the  purpose  of  testing  the  perfonnance  of  the  classification 
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♦ 


For  k  =  1 :  #  of  transmissions 


Compute  the  PWE(t)  based  on  the  SNR  criterion  with  the  proper 
probability  update 


Target  realization  of  g  based  on  the  .■"'  hypothesis 


Simulate  measurement  data  yt  =  ^E.Xj.g.  +  Q. 


Figure  6.  The  pseudo-code  block  diagram  of  one  radar  classification  evolution. 


system  by  randomizing  the  key  variables  in  the  experiment.  By  simulating  many 
realizations  of  every  random  variable  in  the  experiment,  the  Monte  Carlo  simulation  will 
have  a  classification  performance  result  that  is  very  close  to  the  true  performance  of  the 
system.  The  pseudo  code  for  the  Monte  Carlo  simulation  is  presented  in  Figure  7.  The 
Monte  Carlo  simulation  computes  the  probability  of  correct  classification  ( Pcc )  for  the  CR 
system  at  a  specific  transmission  energy  level  E s .  For  this  particular  energy  constraint 

level,  the  simulation  first  chooses  one  target  hypothesis  and  generates  X  number  of 

target  realizations  based  on  this  particular  target  hypothesis.  Secondly,  for  each  target 
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realization,  the  system  generates  Y  measurement  noise  realizations.  Third,  for  each  noise 
realization,  the  system  then  performs  one  cycle  of  target  classification  and  records  the  Pcc 
of  the  classification  result.  Next,  the  system  averages  the  Pcc  of  the  classification  results 
for  this  one  hypothesis.  The  system  repeats  the  same  procedure  for  the  rest  of  the  target 
hypotheses  and  averages  the  Pcc  from  all  hypotheses  to  obtain  the  Pcc  for  this  one 
particular  energy  level. 


►  For  kl  =  1:  #  of  energy  levels 


Energy  Level  =  Energy(kl) 


For  k2  •  1:  S  of  hypotheses  (M) 


Error  =  0 


For  k3  =  1:  U  of  target  realizations 


Generate  a  target  based  on  hypoihesisk2 


For  k4  =  1 :  £  of  noise  realizations 


Perform  target  classification  based  on  one  noise  realization 


If  wrong,  increment  Enor  by  1 


End 


Error  Rate  =  Enor  (=  of  target  realizations)  (  =  of  noise  realization) 


End 


End 


Overall  Error  Rate 'S  Energy(kl)  =  Error  Rate  =  of  hypotheses 


End 


Figure  7.  An  example  of  the  Monte  Carlo  simulation  used  in  the  experiment. 
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For  the  overall  experiment,  the  Monte  Carlo  simulation  system  uses  X  =  50  for 
the  number  of  target  realizations  and  7  =  10  for  the  number  of  noise  realizations.  The 
energy  levels  Es  for  this  research  are  defined  on  a  logarithmic  scale  ranged  from  10  2  to 

10  energy  units.  At  the  end,  an  overall  Pcc  vs.  energy  curve  is  plotted  to  evaluate  the 
performance  of  the  CR  system. 

For  comparison  purposes,  this  research  also  uses  a  traditional  wideband  pulse  as 
the  input  illumination  waveform  to  conduct  the  same  experiments  as  the  PWE  waveform 
design  method.  An  overall  Pcc  vs.  energy  curve  is  plotted  for  the  wideband  illumination 
waveform  and  uses  this  result  as  a  benchmark  to  show  the  improvement  of  the  SNR- 
based  PWE  waveform  method. 

The  experimental  setup  for  this  research  was  explained  in  this  chapter.  This 
includes  the  method  for  target  simulation,  matched  waveform  design,  target  classification 
scheme  implementation,  and  the  Monte  Carlo  simulation  setup  for  test  evaluations.  In  the 
next  chapter,  a  representative  example  the  PWE  waveform,  the  simulated  targets  and  the 
results  from  various  test  runs  are  presented  and  explained. 
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VI.  EXPERIMENTAL  RESULTS 


In  order  to  test  the  algorithm  presented  in  this  research,  a  series  of  experiments 
were  carried  out  to  evaluate  its  perfonnance.  In  this  chapter,  first,  a  representative 
example  of  a  PWE  illumination  waveform  is  presented.  Then,  the  test  results  of  the 
experiments  are  presented  along  with  detail  explanations  and  conclusions. 

A.  REPRESENTATIVE  EXAMPLE  OF  A  PWE  ILLUMINATION 
WAVEFORM 

For  the  purpose  of  illustrating  the  PWE  illumination  process,  a  representative 
example  is  presented  in  this  section.  First,  the  proposed  PSVs  of  the  four  target 
hypotheses  are  created  using  the  parameters  listed  in  Table  1.  The  proposed  PDFs  for  the 
target  to  draw  samples  from  are  exponential  PDFs  with  ju= 2 .  The  parameters  for  the 
PSVs  are  ju  and  a ,  where  //  defines  the  center  frequencies  of  the  energy  bands  of  the 
targets  and  a  defines  the  relative  amplitudes  and  bandwidths  of  the  band-pass  waveform 
spectra.  The  number  of  data  samples  used  in  this  test  is  set  to  64. 


Table  1 .  Key  parameters  for  the  target  hypotheses  in  the  example. 


Target  Hypotheses 

Target  PDFs  (PDF  parameters) 

M 

a 

Target  #  1 

Complex  Exponential  (//  =  2) 

0.15 

0.01 

Target  #  2 

Complex  Exponential  ( //  =  2) 

0.2 

0.015 

Target  #  3 

Complex  Exponential  (//  =  2) 

0.25 

0.02 

Target  #  4 

Complex  Exponential  (//  =  2) 

0.3 

0.025 

39 


The  PSVs  of  the  four  target  hypotheses  are  illustrated  in  Figure  8.  By  putting  the 
four  PSVs  on  one  plot,  we  can  see  that  none  of  the  target  spectra  completely  overlap  any 
of  the  other  target  spectra.  All  target  hypotheses  have  different  amplitudes  and 
bandwidths  in  their  PSVs. 


Figure  8.  PSVs  corresponding  to  the  four  target  hypotheses. 

Next,  a  particular  target  realization  from  target  hypothesis  H{  is  generated.  The 

PSV  of  this  target  realization  is  illustrated  in  Figure  9.  Because  the  experiment  only  sets 
the  number  of  samples  to  64,  the  PSV  of  any  stochastic  target  realization  is  not  very 
smooth  in  nature.  However,  the  generated  target  still  possesses  the  general  properties  of 
the  target  hypothesis  in  terms  of  total  energy  and  band-pass  frequencies.  The  target  in 
Figure  9  and  target  hypothesis  Hx  in  Figure  8  both  have  center  frequency  0.2  unit  on  a 
nonnalized  frequency  scale.  The  total  energy  of  both  PSVs  is  equal  to  1  energy  unit. 
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Figure  9.  PSV  of  one  target  realization  (Target  Hypothesis  #  1). 


The  PWE  illumination  waveform  is  created  by  the  method  described  in  the 
experimental  setup.  The  PSVs  of  the  first  and  the  tenth  illumination  waveforms  for  both 
the  PWE  illumination  waveform  and  the  wideband  illumination  waveforms  are  illustrated 
in  Figure  10.  From  the  PSV  plots  we  can  see  that  for  the  initial  transmission,  the  PWE 
illumination  wavefonn  possesses  the  properties  of  all  four  target  hypotheses.  Then,  at  the 
tenth  transmission,  the  PWE  waveform  converges  and  focuses  all  of  its  energy  towards 
the  frequency  bands  in  which  target  hypothesis  H]  resides.  The  wideband  method  stays 

the  same  for  all  ten  transmissions.  From  this  example,  we  can  see  that  using  the  PWE 
illumination  waveforms  results  in  a  very  efficient  distribution  of  the  limited  illumination 
energy  toward  the  high  energy  parts  of  the  PSV.  When  a  wideband  pulse  is  used,  most 
energy  is  wasted  on  the  frequency  bands  where  the  target  spectra  have  very  little  or  zero 
energy. 
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Figure  10.  The  initial  and  tenth  transmission  snapshots  of  the  PSVs  for  PWE 
illumination  waveform  and  the  wideband  illumination  waveform:  (a)  PSV  of  the  initial 
PWE  illumination  waveform;  (b)  PSV  of  the  initial  wideband  illumination  waveform;  (c) 
PSV  of  the  tenth  PWE  illumination  waveform;  (d)  PSV  of  the  tenth  wideband 

illumination  waveform. 


The  time  waveforms  of  the  PWE  illumination  waveform  at  the  tenth  transmission 
are  illustrated  in  Figure  11.  From  the  time  waveforms  of  the  real  and  imaginary  values  in 
the  PWE  illumination  waveform,  we  can  see  a  rough  looking  weighted  sum  of  sinusoidal 
waveforms.  This  corresponds  to  the  fact  that  the  PWE  waveforms  are  generated  based  on 
an  eigenfunction  solution,  which  produces  sinusoidal  functions. 
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Figure  1 1 .  Time  waveform  of  the  PWE  illumination  wavefonn  at  the  10th 

transmission.  Top  and  bottom  are  the  real  values  and  imaginary  values  of  the  PWE 

illumination  waveform,  respectively. 


This  concludes  the  demonstration  of  the  PWE  illumination  waveform  process. 
Next,  the  experimental  results  from  various  tests  are  presented  and  explained. 

B.  CLASSIFICATION  RESULTS  FOR  FOUR  GAUSSIAN  TARGET 
CLASSES 

The  first  experimental  test  set  consists  of  four  Gaussian  distributed  complex 
targets  with  different  PSVs.  This  test  is  conducted  mainly  to  evaluate  the  performance  of 
the  new  classification  algorithm  on  Gaussian  target  classes.  The  four  proposed  PSVs  are 
illustrated  in  Figure  12.  Each  PSV  is  a  mixture  of  two  Gaussian  distribution  functions. 
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None  of  the  target  PSVs  completely  overlaps  other  target  PSVs.  The  number  of  data 
samples  used  in  this  test  is  set  to  64.  The  specific  parameters  of  the  four  different  target 
hypotheses  are  listed  in  Table  2. 


Table  2.  Parameter  setup  of  the  four  target  hypotheses  for  the  first  experiment. 


Target  Hypotheses 

Target  PDFs  (PDF  parameters) 

M 

<j 

Target  #  1 

Complex  Gaussian)//  =  0,cr  =  1) 

0.15 

0.01 

Target  #  2 

Complex  Gaussian)//  =  0,er  =  1) 

0.2 

0.015 

Target  #  3 

Complex  Gaussian)//  =  0,er  =  1) 

0.25 

0.02 

Target  #  4 

Complex  Gaussian)//  =  0,cr  =  1) 

0.3 

0.025 

The  parameter  //  defines  the  band-pass  center  frequencies  of  the  PSVs,  and  the 
parameter  a  defines  the  relative  amplitudes  and  bandwidths  of  the  energy  bands  of  the 
PSVs. 

An  example  of  one  set  of  target  realizations  based  on  the  proposed  PDFs  and 
PSVs  is  illustrated  in  Figure  13.  Because  the  experiment  only  has  the  number  of  samples 
set  to  64,  the  spectrum  of  any  stochastic  target  realization  is  not  smooth  in  nature. 
However,  the  generated  target  still  possesses  the  general  properties  of  the  target 
hypothesis  proposed  in  Figure  12  in  tenns  of  total  energy  and  band-pass  frequencies.  The 
spectra  presented  in  Figure  13  conform  to  the  respective  PSVs  of  the  target  hypotheses. 

The  PDF  estimates  of  one  set  of  target  realizations  based  on  the  proposed  PDFs 
and  PSVs  are  illustrated  in  Figure  14.  The  PDF  estimates  are  computed  using  a  univariate 
Kernel  Density  Estimator  (KDE)  as  described  in  [19].  The  PDFs  of  both  the  real  and 
imaginary  values  in  all  four  targets  have  a  rough  looking  Gaussian  shape  with  zero  mean 
and  unity  variance,  which  conform  to  the  proposed  Gaussian  PDF. 
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Figure  12.  An  example  of  the  proposed  PSVs  of  four  target  hypotheses:  (a)  target 
hypothesis  one;  (b)  target  hypothesis  two;  (c)  target  hypothesis  three;  (d)  target 

hypothesis  four. 


(a) 


(b) 


(c) 


(d) 


Figure  13.  The  spectra  of  one  realization  of  targets  from  the  first  target  hypothesis  set 
based  on  four  proposed  PSVs  and  a  Gaussian  PDF:  (a)  target  hypothesis  one;  (b)  target 
hypothesis  two;  (c)  target  hypothesis  three;  (d)  target  hypothesis  four. 
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Figure  14.  PDF  estimates  of  the  four  target  realizations.  The  PDF  estimates  of  the 
real  values  of  the  four  target  realizations  are  presented  on  the  left  column  and  the  PDF 
estimates  of  the  imaginary  values  of  the  four  target  realizations  are  presented  on  the  right 

column. 


The  time  waveforms  of  one  set  of  target  realizations  are  illustrated  in  Figure  15. 
From  the  time  waveforms,  we  can  see  the  correlations  in  the  data  samples,  which  also 
indicate  that  the  data  samples  in  these  targets  are  not  independent  due  to  the  ZMNL 
transformation  process. 

The  normal  distribution  plots  of  the  first  two  targets  from  one  set  of  target 
realizations  are  presented  in  Figure  16.  The  nonnal  distribution  plot  is  used  for  evaluating 
the  Gaussianity  of  a  signal.  If  the  data  samples  in  a  signal  are  Gaussian  distributed,  then 
the  data  samples  align  to  the  red  diagonal  line  in  the  normal  distribution  plot.  From  these 
nonnal  distribution  plots,  we  can  see  that  the  data  samples  from  first  two  target 
realizations  are  approximately  Gaussian  distributed. 
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Figure  15.  Time  waveform  of  the  four  target  hypotheses.  Left  column:  real  values  of 
target  realizations;  right  column:  imaginary  values  of  the  target  realizations. 
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Figure  16.  Normal  distribution  plots  of  the  first  two  targets  in  the  target  realization 
set.  Left  column:  normal  distribution  plots  of  the  real  values  of  the  first  two  targets;  right 
column:  normal  distribution  plots  of  the  imaginary  values  of  the  first  two  targets. 
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The  normal  distribution  plots  of  the  last  two  targets  from  one  set  of  target 
realizations  are  presented  in  Figure  17.  From  these  normal  distribution  plots,  we  see  that 
the  data  samples  from  last  two  target  realizations  are  approximately  Gaussian  distributed. 


Re[g3] 


Figure  17.  Normal  distribution  plots  of  the  first  two  targets  in  the  target  realization 
set.  Left  column:  normal  distribution  plots  of  the  real  values  of  the  first  two  targets;  right 
column:  normal  distribution  plots  of  the  imaginary  values  of  the  first  two  targets. 


The  Monte  Carlo  experimental  result  for  one  radar  transmission  is  presented  in 
Figure  18.  The  Monte  Carlo  simulation  is  set  to  generate  50  target  realizations  and  ten 
noise  realizations  for  each  target  realization.  The  number  of  the  online  target  realizations 
for  the  expected  value  estimation  is  set  to  Ng  =  40 .  For  four  target  hypotheses  that  are 

Gaussian  distributed,  the  PWE  illumination  wavefonn  outperforms  the  traditional 
wideband  illumination  waveform  and  saves  approximately  6  dB  in  transmission  energy 
when  the  Pcc  is  90%.  Another  way  to  view  the  performance  gain  is  that  for  a  given  mid- 
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range  transmission  energy,  we  achieve  approximately  a  100%  increase  in  Pcc.  At  very 
low  energy  levels,  the  receiver  can  no  longer  obtain  any  useful  information  from  the 
received  data  since  the  AWGN  in  the  receiver  becomes  the  dominant  component.  At  very 
high  energy  levels,  any  waveform  is  a  good  waveform;  thus,  the  wideband  illumination 
waveform  is  just  as  good  as  the  PWE  illumination  waveform. 


Figure  18.  Classification  performance  curves  of  the  optimized  SNR  illumination 
wavefonn  and  a  traditional  wideband  waveform  with  the  implementation  of  the  NGCCR 
system  for  one  transmission.  Target  distributions:  Gaussian  x  4.  Monte  Carlo  setup:  50 
target  realizations;  ten  noise  realizations.  NGCCR  system  setup:  40  online  target 
realizations  for  the  ensemble  averaging. 


The  Monte  Carlo  experimental  result  for  ten  transmissions  is  presented  in  Figure 
19.  The  Monte  Carlo  simulation  is  set  to  generate  50  target  realizations  and  ten  noise 
realizations  for  each  target  realization.  The  number  of  the  online  target  realizations  for 
the  expected  value  estimation  is  set  to  Ng  =  40  .  The  Pcc  for  both  illumination  waveforms 

have  improved  in  this  case.  For  four  target  hypotheses  that  are  Gaussian  distributed,  the 
PWE  illumination  wavefonn  outperforms  the  traditional  wideband  illumination 
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waveform  and  saves  approximately  5  dB  in  transmission  energy  when  the  Pcc  is  90%. 
Another  way  to  view  the  performance  gain  is  that  for  a  given  mid-range  transmission 
energy,  we  achieve  approximately  a  100%  increase  in  Pcc.  At  very  low  energy  levels,  the 
receiver  can  no  longer  obtain  any  useful  information  from  the  received  data  since  the 
AWGN  in  the  receiver  is  the  dominant  component.  At  very  high  energy  levels,  the 
wideband  illumination  waveform  has  a  100%  Pcc.  However,  the  PWE  illumination 
waveform  only  has  a  near  100%  Pcc.  One  of  the  reasons  may  be  that  more  Monte  Carlo 
averaging  is  needed.  For  the  majority  of  the  time,  the  PWE  illumination  wavefonn  is  still 
the  better  choice. 


Figure  19.  Classification  performance  curves  of  the  optimized  SNR  illumination 
wavefonn  and  a  traditional  wideband  waveform  with  the  implementation  of  the  NGCCR 
system  for  ten  transmissions.  Target  distributions:  Gaussian  x  4.  Monte  Carlo  setup:  50 
target  realizations;  ten  noise  realizations.  NGCCR  system  setup:  40  online  target 
realizations  for  the  ensemble  averaging. 
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c. 


CLASSIFICATION  RESULTS  FOR  FOUR  EXPONENTIAL  TARGET 
CLASSES 


The  second  experimental  test  set  consists  of  four  exponentially  distributed 
complex  targets  with  different  PSVs.  This  test  is  conducted  mainly  to  evaluate  the 
performance  of  the  new  classification  algorithm  on  non-Gaussian  target  classes.  The  four 
proposed  PSVs  are  illustrated  in  Figure  20.  Each  PSV  is  a  mixture  of  two  Gaussian 
distribution  functions.  None  of  the  target  PSVs  completely  overlaps  other  target  PSVs. 
The  number  of  data  samples  used  in  this  test  is  set  to  64.  The  specific  parameters  of  the 
four  different  target  hypotheses  are  listed  in  Table  3. 


Table  3.  Parameter  setup  of  the  four  target  hypotheses  for  the  second  experiment. 


Target  Hypotheses 

Target  PDFs  (PDF  parameters) 

M 

a 

Target  #  1 

Complex  Exponential  (//  =  2) 

0.15 

0.01 

Target  #  2 

Complex  Exponential  (//  =  2) 

0.2 

0.015 

Target  #  3 

Complex  Exponential  ( //  =  2) 

0.25 

0.02 

Target  #  4 

Complex  Exponential  (//  =  2) 

0.3 

0.025 

The  parameter  ju  defines  the  band-pass  center  frequencies  of  the  PSVs,  and  the 
parameter  a  defines  the  relative  amplitudes  and  bandwidths  of  the  energy  bands  of  the 
PSVs. 

An  example  of  one  set  of  target  realizations  based  on  the  proposed  PDFs  and 
PSVs  is  illustrated  in  Figure  21.  Because  the  experiment  only  has  64  samples,  the 
spectrum  of  any  stochastic  target  realization  is  not  very  smooth  in  nature.  However,  the 
generated  target  still  possesses  the  general  properties  of  the  target  hypothesis  proposed  in 
Figure  20  in  terms  of  total  energy  and  band-pass  frequencies.  The  spectra  presented  in 
Figure  21  conform  that  the  targets  possess  the  proposed  PSVs  of  the  target  hypotheses. 
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Figure  20.  The  four  proposed  PSVs  for  the  target  hypotheses  in  the  second 
experimental  test  set:  (a)  target  hypothesis  one;  (b)  target  hypothesis  two;  (c)  target 
hypothesis  three;  (d)  target  hypothesis  four. 


(a) 


(b) 


(c) 


(d) 


Figure  2 1 .  The  spectra  of  one  realization  of  targets  from  the  first  target  hypothesis  set 
based  on  four  proposed  PSVs  and  an  Exponential  PDF:  (a)  target  hypothesis  one;  (b) 
target  hypothesis  two;  (c)  target  hypothesis  three;  (d)  target  hypothesis  four. 

52 


The  PDF  estimates  of  one  set  of  target  realizations  based  on  the  proposed  PDFs 
and  PSVs  are  illustrated  in  Figure  22.  The  PDFs  of  both  the  real  and  imaginary  values  in 
all  four  target  realizations  have  a  rough  looking  exponential  distribution  with  parameter 
ju  =  2  ,  which  conforms  to  the  proposed  exponential  PDF. 
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Figure  22.  PDF  estimates  of  one  particular  realization  of  the  second  target  hypothesis 
set.  The  four  targets  are  exponentially  distributed.  The  real  values  of  the  targets  are 
presented  on  the  left  column  and  the  imaginary  values  of  the  targets  are  presented  on  the 

right  column. 


The  time  waveforms  of  one  set  of  target  realizations  are  illustrated  in  Figure  23. 
From  the  time  waveforms,  we  can  see  the  correlation  in  the  data  samples,  which  also 
indicate  that  the  data  samples  in  these  targets  are  not  independent  due  to  the  ZMNL 
transformation  process. 
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Figure  23.  Time  waveforms  of  one  particular  realization  of  the  second  target 
hypothesis  set.  The  real  values  of  the  targets  are  presented  on  the  left  column  and  the 
imaginary  values  of  the  targets  are  presented  on  the  right  column. 


The  normal  distribution  plots  of  the  first  two  targets  from  one  set  of  target 
realizations  are  presented  in  Figure  24.  From  these  normal  distribution  plots,  we  can  see 
that  the  data  samples  do  not  align  with  the  red  diagonal  line  in  the  normal  distribution 
plot,  which  implies  that  first  two  targets  realizations  are  not  Gaussian  distributed. 

The  normal  distribution  plots  of  the  last  two  targets  from  one  set  of  target 
realizations  are  presented  in  Figure  25.  From  these  normal  distribution  plots,  we  can  see 
that  the  data  samples  do  not  align  with  the  red  diagonal  line  in  the  normal  distribution 
plot,  which  conforms  that  last  two  target  realizations  are  not  Gaussian  distributed. 
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Figure  24.  Normal  distribution  plot  of  one  particular  realization  of  the  second  target 
hypothesis  set.  The  plots  for  the  real  values  of  the  first  two  targets  are  presented  on  the 
left  column  and  the  plots  for  the  imaginary  values  of  the  first  two  targets  are  presented  on 

the  right  column. 


The  Monte  Carlo  experimental  result  for  one  radar  transmission  is  presented  in 
Figure  26.  The  Monte  Carlo  simulation  is  set  to  generate  50  target  realizations  and  ten 
noise  realizations  for  each  target  realization.  The  number  of  the  online  target  realizations 
for  the  expected  value  estimation  is  set  to  Ng  =  40 .  For  four  target  hypotheses  that  are 

exponentially  distributed,  the  PWE  illumination  waveform  outperforms  the  traditional 
wideband  illumination  waveform  by  approximately  10  dB  in  transmission  energy  when 
Pcc  is  90%.  Another  way  to  view  the  performance  gain  is  that  for  a  given  mid-range 
transmission  energy,  we  achieve  approximately  a  100%  increase  in  Pcc.  At  very  low 
energy  levels,  the  receiver  can  no  longer  obtain  any  useful  information  from  the  received 
data  since  the  AWGN  in  the  receiver  is  the  dominant  component.  At  very  high  energy 
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levels,  any  waveform  is  a  good  waveform;  thus,  the  wideband  illumination  wavefonn  is 
just  as  good  as  the  PWE  illumination  waveform  at  high  energy  levels. 


Re[g3] 


imlg-sl 


Figure  25.  Normal  distribution  plot  of  one  particular  realization  of  the  second  target 
hypothesis  set.  The  plots  for  the  real  values  of  the  last  two  targets  are  presented  on  the 
left  column  and  the  plots  for  the  imaginary  values  of  the  last  two  targets  are  presented  on 

the  right  column. 


The  Monte  Carlo  experimental  result  for  ten  radar  transmissions  is  presented  in 
Figure  27.  The  Monte  Carlo  simulation  is  set  to  generate  50  target  realizations  and  ten 
noise  realizations  for  each  target  realization.  The  number  of  the  online  target  realizations 
for  the  expected  value  estimation  is  set  to  Ng  =  40  .  The  Pcc  for  both  illumination 

waveforms  have  improved  in  this  case.  For  four  target  hypotheses  that  are  exponentially 
distributed,  the  PWE  illumination  waveform  outperforms  the  traditional  wideband 
illumination  waveform  by  approximately  4  dB  in  transmission  energy  when  Pcc  is  90%. 
Another  way  to  view  the  performance  gain  is  that  for  a  given  mid-range  transmission 
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energy,  we  achieve  approximately  a  100%  increase  in  Pcc.  At  very  low  energy  levels,  the 
receiver  can  no  longer  obtain  any  useful  information  from  the  received  data  since  the 
AWGN  in  the  receiver  is  the  dominant  component.  At  very  high  energy  levels,  the 
wideband  illumination  waveform  has  a  100%  Pcc ■  However,  the  PWE  illumination 
waveform  only  has  a  near  100%  Pcc.  One  of  the  reasons  may  be  that  more  Monte  Carlo 
averaging  is  needed.  For  the  majority  of  the  time,  the  PWE  illumination  waveform  is  still 
the  better  choice.  Also,  the  Pcc  curves  for  the  non-Gaussian  distributed  targets  are  very 
similar  to  the  Pcc  curves  for  the  Gaussian  distributed  targets.  This  result  shows  the 
effectiveness  of  the  NGCCR  detection  algorithm  for  classifying  non-Gaussian  targets. 


Figure  26.  Classification  performance  curves  of  the  optimized  SNR  illumination 
waveform  and  a  traditional  wideband  waveform  with  the  implementation  of  the  NGCCR 
system  for  one  transmission.  Target  distributions:  exponential  x  4.  Monte  Carlo  setup:  50 
target  realizations;  ten  noise  realizations.  NGCCR  system  setup:  40  online  target 
realizations  for  the  ensemble  averaging. 
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Figure  27.  Classification  performance  curves  of  the  optimized  SNR  illumination 
waveform  and  a  traditional  wideband  waveform  with  the  implementation  of  the  NGCCR 
system  for  ten  transmissions.  Target  distributions:  exponential  x  4.  Monte  Carlo  setup:  50 
target  realizations;  ten  noise  realizations.  NGCCR  system  setup:  40  online  target 
realizations  for  the  ensemble  averaging. 


D.  CLASSIFICATION  RESULTS  FOR  FOUR  NON-GAUSSIAN  TARGET 

CLASSES 

The  last  experimental  test  set  consists  of  four  non-Gaussian  distributed  complex 
targets  with  different  PSVs.  The  non-Gaussian  PDFs  are  Rayleigh,  exponential,  Gamma, 
and  log-normal  distributions.  This  test  is  conducted  mainly  to  evaluate  the  performance 
of  the  new  algorithm  for  non-Gaussian  target  classes  with  four  different  distributions. 
The  four  proposed  PSVs  are  illustrated  in  Figure  28.  Each  PSV  is  a  mixture  of  two 
Gaussian  distribution  functions.  None  of  the  target  PSVs  completely  overlaps  other  target 
PSVs.  The  number  of  data  samples  used  in  this  test  is  set  to  64.  The  specific  parameters 
of  the  four  different  target  hypotheses  are  listed  in  Table  4. 
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Table  4.  Parameter  setup  of  the  four  target  hypotheses  for  the  third  experiment. 


Target  Hypotheses 

Target  PDFs  (PDF  parameters) 

M 

<j 

Target  #  1 

Complex  Rayleigh  (cr  =  10) 

0.15 

0.01 

Target  #  2 

Complex  Exponential  (//  =  2) 

0.2 

0.015 

Target  #  3 

Complex  Gamma  (k  =  2, 6  =  2) 

0.25 

0.02 

Target  #  4 

Complex  Log-Normal  (//  =  0,  a  -  1) 

0.3 

0.025 

Figure  28.  The  four  proposed  PSVs  for  the  target  hypotheses  in  the  last  experimental 
test  set:  (a)  target  hypothesis  one;  (b)  target  hypothesis  two;  (c)  target  hypothesis  three; 

(d)  target  hypothesis  four. 
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The  parameter  //  defines  the  band-pass  center  frequencies  of  the  PS  Vs,  and  the 
parameter  a  defines  the  relative  amplitudes  and  bandwidths  of  the  energy  bands  of  the 
PSVs. 

An  example  of  one  set  of  target  realizations  based  on  the  proposed  PDFs  and 
PSVs  is  illustrated  in  Figure  29.  Because  the  experiment  only  has  64  samples,  the 
spectrum  of  any  stochastic  target  realization  is  not  smooth  in  nature.  However,  the 
generated  target  still  possesses  the  general  properties  of  the  target  hypothesis  proposed  in 
Figure  28  in  terms  of  total  energy  and  band-pass  frequencies.  The  spectra  presented  in 
Figure  29  conform  with  the  specified  PSVs  of  the  target  hypotheses. 


(a) 


(b) 


(c) 


(d) 


Figure  29.  The  spectra  of  one  realization  of  targets  from  the  first  target  hypothesis  set 
based  on  four  proposed  PSVs  and  four  different  non-Gaussian  PDFs:  (a)  target 
hypothesis  one;  (b)  target  hypothesis  two;  (c)  target  hypothesis  three;  (d)  target 

hypothesis  four. 


The  PDF  estimates  of  one  set  of  target  realizations  based  on  the  proposed  PDFs 
and  PSVs  are  illustrated  in  Figure  30.  The  PDFs  of  both  the  real  and  imaginary  values  in 

the  first  target  have  a  rough  looking  Rayleigh  distribution  with  parameter  cr- 10;  the 
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PDFs  of  both  the  real  and  imaginary  values  in  the  second  target  have  a  rough  looking 
Exponential  distribution  with  parameter  /u  =  2  ;  the  PDFs  of  both  the  real  and  imaginary 

values  in  the  third  target  have  a  rough  looking  Gamma  distribution  with  parameter  k  =  2 
and  0  =  2;  the  PDFs  of  both  the  real  and  imaginary  values  in  the  fourth  target  have  a 
rough  looking  Log-Normal  distribution  with  parameter  /u  =  0  and  a  =  1 .  These  PDF 
estimates  show  that  these  target  realizations  all  approximate  the  proposed  non-Gaussian 
PDFs. 
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Figure  30.  PDF  estimates  of  one  particular  realization  of  the  last  target  hypothesis  set. 

The  proposed  PDFs  for  the  four  targets  are  Rayleigh,  exponential,  Gamma,  and  log¬ 
normal  PDFs  respectively.  The  real  values  of  the  targets  are  presented  on  the  left  column 
and  the  imaginary  values  of  the  targets  are  presented  on  the  right  column. 
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The  time  waveforms  of  one  set  of  target  realizations  are  illustrated  in  Figure  31. 
From  the  time  waveforms,  we  can  see  the  correlation  in  the  data  samples,  which  also 
indicate  that  the  data  samples  in  these  targets  are  not  independent  due  to  the  ZMNL 
transformation  process. 
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Figure  3 1 .  Time  waveforms  of  one  particular  realization  of  the  last  target  hypothesis 
set.  The  real  values  of  the  targets  are  presented  on  the  left  column  and  the  imaginary 
values  of  the  targets  are  presented  on  the  right  column. 


The  normal  distribution  plots  of  the  first  two  targets  from  one  set  of  target 
realizations  are  presented  in  Figure  32.  From  these  normal  distribution  plots,  we  can  see 
that  the  data  samples  do  not  align  to  the  red  diagonal  line  in  the  normal  distribution  plot, 
which  confirms  that  first  two  target  realizations  are  not  Gaussian  distributed. 

The  normal  distribution  plots  of  the  last  two  targets  from  one  set  of  target 
realizations  are  presented  in  Figure  33.  From  these  normal  distribution  plots,  we  can  see 
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that  the  data  samples  do  not  align  to  the  red  diagonal  line  in  the  normal  distribution  plot, 
which  confirms  that  last  two  target  realizations  are  not  Gaussian  distributed. 


Figure  32.  Normal  distribution  plot  of  one  particular  realization  of  the  last  target 
hypothesis  set.  The  plots  for  the  real  values  of  the  first  two  targets  are  presented  on  the 
left  column  and  the  plots  for  the  imaginary  values  of  the  first  two  targets  are  presented  on 

the  right  column. 


The  Monte  Carlo  experimental  result  for  one  radar  transmission  is  presented  in 
Figure  34.  The  Monte  Carlo  simulation  is  set  to  generate  50  target  realizations  and  ten 
noise  realizations  for  each  target  realization.  The  number  of  the  online  target  realizations 
for  the  expected  value  estimation  is  set  to  N  =  40 .  For  four  target  hypotheses  that  are 

Rayleigh,  exponential,  Gamma,  and  log-normal  distributed,  respectively,  the  PWE 
illumination  waveform  outperforms  the  traditional  wideband  illumination  waveform  by 
approximately  6  dB  in  transmission  energy  when  Pcc  is  90%.  Another  way  to  view  the 
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performance  gain  is  that  for  a  given  mid-range  transmission  energy,  we  achieve 
approximately  a  100%  increase  in  Pcc.  At  very  low  energy  levels,  the  receiver  can  no 
longer  obtain  any  useful  information  from  the  received  data  since  the  AWGN  in  the 
receiver  is  the  dominant  component.  At  very  high  energy  levels,  any  waveform  is  a  good 
waveform;  thus,  the  wideband  illumination  waveform  is  just  as  good  as  the  PWE 
illumination  waveform  at  high  energy  levels.  Note  that  for  this  set  of  non-Gaussian 
targets,  it  appears  that  Pcc  for  the  PWE  waveform  does  not  reach  the  100%  level  after 
only  one  transmission  (at  least  not  over  the  range  of  energy  levels  plotted).  However, 
after  ten  iterations,  we  see  that  the  Pcc  for  the  PWE  waveform  reaches  100%  (see  Figure 
35). 


Figure  33.  Normal  distribution  plot  of  one  particular  realization  of  the  last  target 
hypothesis  set.  The  plots  for  the  real  values  of  the  last  two  targets  are  presented  on  the 
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left  column  and  the  plots  for  the  imaginary  values  of  the  last  two  targets  are  presented  on 

the  right  column. 


Figure  34.  Classification  performance  curves  of  the  optimized  SNR  illumination 
waveform  and  a  traditional  wideband  waveform  with  the  implementation  of  the  NGCCR 
system  for  one  transmission.  Target  distributions:  Rayleigh,  exponential,  Gamma,  log¬ 
normal.  Monte  Carlo  setup:  50  target  realizations;  ten  noise  realizations.  NGCCR  setup: 
40  online  target  realizations  for  the  ensemble  averaging. 


The  Monte  Carlo  experimental  result  for  ten  radar  transmissions  is  presented  in 
Figure  35.  The  Monte  Carlo  simulation  is  set  to  generate  50  target  realizations  and  ten 
noise  realizations  for  each  target  realization.  The  number  of  the  online  target  realizations 
for  the  expected  value  estimation  is  set  to  Ng  =40  .  The  Pcc  for  both  illumination 

waveforms  have  improved  in  this  case.  For  four  target  hypotheses  that  are  Rayleigh, 
exponential,  Gamma,  and  log-normal  distributed,  respectively,  the  PWE  illumination 
waveform  outperforms  the  traditional  wideband  illumination  waveform  by  approximately 
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4  dB  in  transmission  energy  when  Pcc  is  90%.  Another  way  to  view  the  performance  gain 
is  that  for  a  given  mid-range  transmission  energy,  we  achieve  approximately  a  100% 
increase  in  Pcc.  At  very  low  energy  levels,  the  receiver  can  no  longer  obtain  any  useful 
information  from  the  received  data  since  the  AWGN  in  the  receiver  is  the  dominant 
component.  At  very  high  energy  levels,  the  wideband  illumination  waveform  has  a  100% 
Pcc .  However,  the  PWE  illumination  waveform  only  has  a  near  100%  Pcc ■  One  of  the 
reasons  may  be  that  more  Monte  Carlo  averaging  is  needed.  For  majority  of  the  time,  the 
PWE  illumination  waveform  is  still  the  better  choice.  Again,  the  Pcc  curves  for  the  non- 
Gaussian  distributed  targets  are  very  similar  to  the  Pcc  curves  for  the  Gaussian  distributed 
targets.  This  result  also  shows  the  effectiveness  of  the  NGCCR  detection  algorithm  for 
classifying  non-Gaussian  targets. 


Figure  35.  Classification  performance  curves  of  the  optimized  SNR  illumination 
waveform  and  a  traditional  wideband  waveform  with  the  implementation  of  the  NGCCR 
system  for  ten  transmissions.  Target  distributions:  Rayleigh,  exponential,  Gamma,  log¬ 
normal.  Monte  Carlo  setup:  50  target  realizations;  ten  noise  realizations.  NGCCR  setup: 
40  online  target  realizations  for  the  ensemble  averaging. 
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Overall,  we  can  draw  conclusion  from  the  results  from  these  experimental  test  runs  that 
the  PWE  illumination  wavefonn  outperforms  the  traditional  wideband  illumination 
waveform  in  an  energy  constraint  environment  due  to  the  efficient  use  of  the  energy  in 
the  illumination  waveforms.  The  new  proposed  NGCCR  target  classification  scheme 
performs  well  with  both  Gaussian  and  non-Gaussian  distributed  complex  stochastic 
targets  compared  to  wideband  illumination  wavefonn.  In  the  next  chapter,  we  present 
conclusions  and  discuss  topics  for  future  research. 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


Romero  et  al.  recently  developed  a  new  multiple  hypothesis  cognitive  radar  target 
classification  algorithm  (dubbed  here  the  GCCR  system)  that  exploits  the  spectral 
sparsity  of  correlated  narrowband  target  responses  to  achieve  significant  performance 
improvements  over  traditional  radars  that  use  wideband  illumination  pulses. 

The  theory  of  the  GCCR  system  to  deal  effectively  with  arbitrary  zero-mean  non- 
Gaussian  distributed  target  responses  was  generalized  in  this  thesis.  The  new  algorithm  is 
dubbed  the  NGCCR  system.  As  discussed  in  Chapter  3,  the  GCCR  system  exploits 
Gaussian-motivated  closed-form  expressions  to  evaluate  the  joint  density  of  the  multiple 
measurements  in  integral  form 


pkt(yi>yi> 
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The  NGCCR  system  relaxes  the  Gaussian  target  assumption.  In  this  case,  we  cannot 
evaluate  the  expected  value  in  integral  form  since  there  exists  no  closed-form  expression 
for  the  integrand.  The  desired  generalization  was  achieved  via  two  key  contributions. 
First,  the  joint  density  of  the  multiple  measurements  is  evaluated  numerically  by 
replacing  the  integral  with  an  ensemble  averaging  operation  to  evaluate  the  expected 
value  of  the  joint  condition  posterior  density  of  the  measurements  conditioned  on  the 
target  realization  given  by 


(57) 


Next,  a  powerful  new  statistical  sampling  algorithm  by  Nichols  et  al.  and  a  kernel  density 
estimator  were  applied  to  draw  complex  target  samples  from  target  distributions  specified 
by  both  desired  power  spectral  densities  and  arbitrary  desired  probability  density 
functions.  These  realizations  are  then  used  to  form  the  ensemble  of  posterior  density 
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values  employed  in  the  ensemble  average  operation.  The  NGCCR  system  requires  prior 
knowledge  of  the  target  classes  in  the  form  of  target  response  measurements  or 
simulations  based  upon  measurements.  Given  these  measurements,  it  is  possible  to 
estimate  the  PDFs  of  the  target  classes  using  an  appropriate  PDF  estimator  (i.e.,  the 
kernel  density  estimator  [19].  For  this  thesis,  these  PDF  estimates  were  computed  to  aid 
visual  analysis  of  the  data  and  to  gain  insight;  however,  they  were  not  used  in  the 
NGCCR  system  because  they  were  not  required  for  the  implementation.  In  an  operational 
scenario  in  which  the  target  distributions  are  not  known  a  priori,  these  PDFs  can  be 
estimated  and  used  as  part  of  the  statistical  sampling  algorithm.  Thus,  the  NGCCR 
system  does  not  require  prior  knowledge  of  the  “true”  target  distributions.  If  we  are  given 
representative  measurements  of  the  desired  target  responses  (as  we  are  in  typical 
operational  scenarios),  we  can  estimate  their  PDFs,  then  draw  target  samples  from  their 
associated  distributions,  and  use  those  target  samples  to  evaluate  the  desired  joint 
densities  required  above.  Once  we  have  these  joint  densities,  we  use  them  to  update  the 
prior  probabilities.  Next,  we  use  the  updated  priors  to  update  the  PWE  illumination  signal, 
which  is  then  transmitted  by  the  radar  to  illuminate  the  targets.  The  user  defines  in 
advance  the  desired  number  of  iterations  to  use.  When  that  number  of  iterations  is 
complete,  the  final  target  classification  decision  is  made. 

Simulations  using  non-Gaussian  targets  demonstrated  very  effective  classification 
algorithm  performance.  As  expected,  this  performance  gain  is  realized  at  the  expense  of 
increased  computational  complexity. 

Recommended  future  research  lies  in  three  main  areas. 

First,  future  research  should  focus  on  finding  ways  to  reduce  the  computational 
complexity  of  the  new  algorithm. 

Second,  in  addition,  our  current  statistical  sampling  algorithm  is  limited  to  being 
able  to  sample  from  distributions  that  are  specified  by  both  a  desired  symmetric  PSV  and 
a  desired  PDF.  Future  work  should  also  focus  on  creating  a  generalized  method  for 
sampling  a  distribution  specified  by  both  a  desired  asymmetric  PSV  and  a  desired  PDF. 
This  will  likely  involve  using  a  phase  retrieval  algorithm  for  evaluating  the  IFT  of  a  PSV. 
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Third,  a  demonstration  of  the  full  end-to-end  process  of  generating  target  samples 
from  a  representative  arbitrarily  distributed  target  process  should  be  conducted.  This 
specific  exercise  was  not  conducted  for  the  thesis.  Assuming  we  are  given  only 
representative  measurements  from  a  target  process,  we  need  to  estimate  the  PDF  of  the 
target  distribution  using  a  kernel  density  estimator.  Next,  we  should  draw  IID  samples 
from  the  target  process  using  a  statistical  sampling  algorithm  such  as  the  Importance 
Sampling  algorithm  or  the  Metropolis-blastings  Sampling  algorithm.  These  IID  samples 
should  next  be  inserted  into  the  Nichols  sampling  algorithm,  along  with  the  desired  PSV. 
The  results  of  the  Nichols  algorithm  will  be  target  samples  that  have  both  the  desired 
arbitrary  PDF  and  the  desired  PSV.  Completion  of  this  exercise  will  demonstrate  that  in 
an  operational  scenario,  no  prior  knowledge  of  the  target  distribution  is  necessary. 
Clearly,  the  gains  made  by  demonstrating  this  capability  will  be  realized  at  the  expense  of 
greatly  increased  computational  complexity. 
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APPENDIX 


This  code  was  written  to  generate  the  four  PS  Vs  of  the  target  hypotheses  and  evaluate  the 
target  generating  function. 


%  Author  Kenan  Wang 

generate  non-Gaussian  colored  spectra  complex  signals  for  the  CR 
experiment 


close  all; 

clear; 

clc; 


64;  %  define  the  number  of  samples  in  the  signal  sequence 


Ps  =  1;  %  define  the  power  of  the  signal 


define  the  frequency  bands  for  the  bandpass  PSD,  the  frequency  is 
normalized  to  0.5  Hz 


mul  =  0.2; 
mu2  =  0.3; 
mu3  =  0.35; 
mu4  =  0.4; 


define  the  bandwidth  of  the  PSD,  different  sigma  will  result  in  different 
amplitude  and  bandwidth  of  the  PSD 


sigmal  =  0.03; 
sigma2  =  0.02; 
sigma3  =  0.025; 
sigma 4  =  0.015; 


define  the  four  pdfs  we  want  the  complex  targets  to  draw  samples  from 
so  far  only  these  four  options  are  allow.  Adding  options  requires  the 
changing  of  code  in  the  function  "generate_targets . m" 
pdf  options  are:  'rayleigh' ,' exponential' ,' gamma' , and  'lognorm' 


pdfl  =  'rayleigh' ; 
pdf 2  =  'exponential'; 
pdf 3  =  'gamma' ; 
pdf 4  =  'lognorm' ; 

define  the  normalized  frequency  domain 
df  =  1/N; 

f  =  -0 . 5 : df : 0 . 5-df  ; 
compute  the  four  proposed  PSDs 

Sxxl  =  (Ps/(2*sqrt(2*pi) *  sigmal ) ) * (exp ( - ( f-mul)  .  A2/ (2  * sigmal A2 ) ) +exp ( - 
( f+mul )  . A2/ (2  * sigmal A2 ) ) ) ; 

Sxx2  =  (Ps/ (2*sqrt(2*pi) *  sigma2 ) ) * (exp ( - ( f-mu2 )  . A2/ (2  *sigma2 A2 ) ) +exp ( - 
(f+mu2) ,A2/ (2*sigma2A2) ) ) ; 
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Sxx3  =  (Ps/ (2*sqrt (2*pi) *sigma3) ) * (exp (- (f-mu3) . A2/ (2*sigma3A2) ) +exp (- 
(f+mu3) .  A2/ (2*sigma3A2) ) ) ; 

Sxx4  =  (Ps/(2*sqrt(2*pi) *  sigma4 ) ) * (exp ( - ( f-mu4 )  . A2/ (2  * sigma 4 A2 ) ) +exp ( - 
( f+mu4 )  . A2/ (2  * sigma 4 A2 ) ) ) ; 


plot  the  four  proposed  PSDs 

figure 

subplot  (221 ) 

plot ( f , Sxxl ) ; 

ylabel ( 'Sxxl (f) ' ) 

xlabel ( 'Normalized  Frequency' ) 

title ( 'Target  Hypothesis  #1') 

subplot  (222 ) 

plot ( f , Sxx2 ) ; 

ylabel (  'Sxx2 (f) ' ) 

xlabel ( 'Normalized  Frequency' ) 

title  ( 'Target  Hypothesis  #2') 

subplot (223 ) 

plot (f, Sxx3) ; 

ylabel  (  'Sxx3 (f) ' ) 

xlabel ( 'Normalized  Frequency' ) 

title  ( 'Target  Hypothesis  #3') 

subplot (224 ) 

plot  ( f , Sxx4 ) ; 

ylabel (  'Sxx4 (f) ' ) 

xlabel ( 'Normalized  Frequency' ) 

title  ( 'Target  Hypothesis  #4') 

simulate  four  target  realizations  based  on  the  proposed  pdf  and  PSDs 


hi  =  generate_targets (  Sxxl, 

%hl  =  hi  -  mean (hi); 

h2  =  generate_targets (  Sxx2, 

%h2  =  h2  -  mean(h2) ; 

h3  =  generate_targets (  Sxx3, 

%h3  =  h3  -  mean (h3) ; 

h4  =  generate_targets (  Sxx4 , 

%h4  =  h4  -  mean(h4); 


pdfl,  N) 
pdf2,  N) 
pdf 3,  N) 
pdf 4,  N) 


+  1 j *generate_targets (  Sxxl, 
+  1 j *generate_targets (  Sxx2, 
+  1 j *generate_targets (  Sxx3, 
+  1 j *generate_targets (  Sxx4  , 


pdfl,  N) 
pdf2,  N) 
pdf 3,  N) 
pdf 4,  N) 


compute  the  PSV  of  the  four  realized  target  and  plot  the  results 

HI  =  abs (fftshift (fft (hi) ) ) . A2; 

H2  =  abs (fftshift  (fft (h2) ))  ,A2; 

H3  =  abs (fftshift (fft (h3) )) ,A2; 

H4  =  abs (fftshift (fft (h4))).A2; 


figure 

subplot  (221 ) 

plot  ( f , HI ) 

ylabel (  'Sxxl (f) ' ) 

xlabel ( 'Normalized  Frequency' ) 

title  ( 'Target  #1') 

subplot (222 ) 

plot ( f , H2 ) 

ylabel (  'Sxx2 (f) ' ) 

xlabel ( 'Normalized  Frequency' ) 

title  ( 'Target  #2') 

subplot (223 ) 

plot (f, H3) 

ylabel ( 'Sxx3 (f) ' ) 

xlabel ( 'Normalized  Frequency' ) 
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title  ( 'Target  #3') 

subplot (224 ) 

plot ( f , H4 ) 

ylabel (  'Sxx4 (f) ' ) 

xlabel ( 'Normalized  Frequency' ) 

title  ( 'Target  #4') 


the  following  code  allow  you  to  see  the  pdf  estimates,  the  time  waveform 
and  the  normplots  of  the  four  complex  target  signals 


[f_hl_real  hil_real] 
[f_hl_imag  hil_imag] 
[f_h2_real  hi2_real] 
[f_h2_imag  hi2_imag] 
[f_h3_real  hi3_real] 
[f_h3_imag  hi3_imag] 
[f_h4_real  hi4_real] 
[f_h4_imag  hi4_imag] 


=  ksdensity (real (hi) ) ; 
=  ksdensity ( imag (hi )) ; 
=  ksdensity (real (h2) ) ; 
=  ksdensity ( imag (h2 )) ; 
=  ksdensity (real (h3) ) ; 
=  ksdensity (imag (h3) ) ; 
=  ksdensity (real (h4) ) ; 
=  ksdensity ( imag (h4 )) ; 


figure 
subplot  (421 ) 

plot (hil_real ,  f_hl_real) 
ylabel ( 'p (real (hi) ) ' ) 
xlabel (  'real (hi ) ' ) 
subplot (422 ) 

plot  (hil_imag,  f_hl_imag) 
ylabel ( 'p (imag (hi ) ) ' ) 
xlabel ( 'imag (hi ) ' ) 
subplot (423) 

plot (hi2_real ,  f_h2_real) 
ylabel (  'p  (real  (h2) ) ' ) 
xlabel (  'real (h2 ) ' ) 
subplot (424 ) 

plot (hi2_imag,  f_h2_imag) 
ylabel ( 'p (imag (h2 ) ) ' ) 
xlabel ( 'imag (h2 ) ' ) 
subplot (425) 

plot (hi3_real ,  f_h3_real) 
ylabel (  'p  (real  (h3) ) ' ) 
xlabel (  'real (h3) ' ) 
subplot (426) 

plot (hi3_imag,  f_h3_imag) 
ylabel ( 'p (imag (h3) ) ' ) 
xlabel ( 'imag (h3) ' ) 
subplot (427 ) 

plot  (hi4_real ,  f_h4_real) 
ylabel (  'p  (real (h4) )  ' ) 
xlabel (  'real (h4 ) ' ) 
subplot (428) 

plot (hi4_imag,  f_h4_imag) 
ylabel ( 'p (imag (h4 ) ) ' ) 
xlabel ( 'imag (h4 ) ' ) 


figure 
subplot  (421 ) 
plot (real (hi) ) 

title ('time  waveform  of  real (hi)') 
ylabel ( 'real (hi ) ' ) 
xlabel (  'n' ) 
subplot  (422 ) 
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plot ( imag (hi ) ) 

title ( 'time  waveform  of  imag (hi)') 

ylabel ( 'imag (hi ) ' ) 

xlabel (  'n' ) 

subplot (423) 

plot  (real  (h2) ) 

title ( 'time  waveform  of  real(h2)') 

ylabel (  'real (h2 ) ' ) 

xlabel  (  'n' ) 

subplot (424 ) 

plot ( imag (h2 ) ) 

title ( 'time  waveform  of  imag(h2)') 

ylabel (  'imag (h2 ) ' ) 

xlabel (  'n' ) 

subplot (425) 

plot (real  (h3) ) 

title ( 'time  waveform  of  real(h3)') 

ylabel (  'real (h3) ' ) 

xlabel ( 'n' ) 

subplot (426) 

plot ( imag (h3) ) 

title ( 'time  waveform  of  imag(h3)') 

ylabel ( 'imag (h3) ' ) 

xlabel (  'n' ) 

subplot (427 ) 

plot (real  (h4) ) 

title ( 'time  waveform  of  real(h4)') 

ylabel ( 'real (h4 ) ' ) 

xlabel ( 'n' ) 

subplot (428) 

plot ( imag (h4 ) ) 

title (  'time  waveform  of  imag(h4)') 
ylabel (  'imag (h4 ) ' ) 
xlabel ( 'n' ) 

figure 
subplot  (221 ) 
normplot (real (hi ) ) 
subplot (222 ) 
normplot ( imag (hi ) ) 
subplot (223 ) 
normplot (real (h2 ) ) 
subplot (224 ) 
normplot ( imag (h2 ) ) 

figure 
subplot  (221 ) 
normplot (real  (h3) ) 
subplot (222 ) 
normplot ( imag (h3) ) 
subplot (223 ) 
normplot (real (h4 ) ) 
subplot (224 ) 
normplot ( imag (h4 ) ) 


save  the  PSV  and  the  names  of  the  pdfs  of  the  targets  for  the  experiment 


Sxx  =  [Sxxl'  Sxx2'  Sxx3'  Sxx4 ' ] ; 
pdf_names  =  {pdfl  pdf2  pdf3  pdf4}; 
save (  'targets .mat' , ' Sxx' , 'pdf_names' ) 
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This  code  was  written  for  the  Monte  Carlo  simulation  of  the  multiple  hypothesis  target 
classification  algorithm.  The  PWE  and  the  wideband  illumination  waveforms  are  used  for 
comparisons. 

%  Author:  Kenan  Wang 

%  Case:  non-Gaussian  stochastic  targets  in  AWGN 
%  Waveform  Design:  Uses  PWE  (probability  weighted  energy) 

%  Target  Recognition  Experiment 
%  Estimate  of  the  maximum  likelihood  Detection 

close  all; 

clear; 

clc; 


ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 


%  INITIALIZATIONS:  The  user  should  choose  parameters  here 


ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 


%  Monte  Carlo  loop  parameters 


target  realization 

=  50 

num  exp 

=  10 

Ntheta 

=  4; 

%  Number  of  target  realizations 
%  Number  of  noise  realizations 
%  No. of  target  hypotheses  (this  is  fixed) 


transmission 
num  ensemble 


=10;  %  Number  of  radar  transmissions 

=40;  %  Number  of  the  target  realizations  used  in  the 

%  online  ensemble  to  calculate  the  expected 
%  value  of  the  difference  between  the  received 
%  data  and  the  target  hypotheses 
%  this  number  is  good  for  10-40  but  above  25 
%  will  see  very  consistent  results 


load  targets; 


%  load  the  PSV  of  the  target  hypotheses 
%  this  .mat  file  contains  the  Sxx  and  the  name 
%  of  the  target  distributions 


Lsignal  =  length (Sxx); 
Ly  =  Lsignal*2-1; 


%  obtain  the  length  of  the  signal 
%  compute  the  length  of  the  convolution  between 
%  the  input  and  the  signal  process 


Pn  =  1;  %  Noise  power,  setting  it  to  1  will  give  TNR  =  lOdB 

var_noise  =  sqrt(Pn/2);  %  compute  the  variance  of  the  noise 

Kn  =  var_noise*eye (Ly ) ;  %  compute  the  covariance  matrix  of  the  noise 


N  =  Lsignal; 


%  rename  the  variable  for  convenience 
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ES  =  energy_levels ( le-2 , 1 0 ) ;  %  define  the  range  of  the  levels  of  energy  for  the 
radar  transmitter 

%  that  we  want  to  run  the  experiment  on 


%  the  fi: 
%  second 

accuracy  =  zeros ( 1 , length (ES )) ; 
vector 

final_result  =  zeros (2, length (ES) ) ; 
result  vector 

dec  =  0 ; 

loop_track  =  0; 

for  the  experiment 


st  argument  is  the  lower  bound  and  the 
argument  is  the  upperbound 

%  Memory  allocation  for  the  accuracy 

%  Memory  allocation  for  the  final 

%  initialize  the  decision  variable 
%  a  variable  to  track  the  num.  of  loops 


%  SNR  CASE  WITH  NO  CLUTTER: 

%  Compute  the  Optimal  Input  Signal  xi (t)  for  each  of  the  Target 
%  Hyptheses  i  =  1,  2,  .  .  .  ,  M,  where  M  =  Nthetha  in  the  code 

%  the  second  argument  indicates  whether  the  the  target  hypotheses  are  in  the 
frequency 

%  domain  or  the  time  domain.  0  for  frequency  domain,  1  for  time  domain 


x_opt  =  opt_wavef orms (Sxx' , 0) ; 


%  Create  a  complex  wideband  input 


delta  =  zeros (Lsignal, 1) ; 

delta  (1)  =  1; 

delta  =  delta  +  lj*delta; 

delta  =  delta/sqrt (delta' *delta) ; 


for  kl  =  1:2  %  perform  this  process  twice,  first  time  is  the  PWE  input  method, 
the  second  time  is  the  wideband  method 

if  (kl  ==  1) 

wavef orm_method  =  'PWE' ;  %  Define  the  waveform  method  used 

else 

wavef orm_method  =  'wideband' ; 

end 


for  k2  =  1: length  (ES) 
energy  level 


%  perform  target  classification  at  each 


Es  =  ES (k2)  ; 

for  k3  =  l:Ntheta 
hypotheses 

err  =  0; 
hypo  =  k3 ; 

hypothesis 


%  Energy  per  transmission 
%  Monte  Carlo  loop  for  the  four  target 

%  initialize  the  error 
%  set  the  hypo  equal  to  the  defined 
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for  k4  =  1 : target_realization  %  Monte  Carlo  loop  for  number  of 
target  realizations 


%  an  algorithm  to  generate  the  non-Gaussian  targets  based 
%  on  the  current  hypothesis 


switch  hypo 
case  1 

h  =  generate_targets (  Sxx(:,l)', 
N)  +  1 j *generate_targets (  Sxx(:,l)',  char (pdf_names ( 1 )  )  ,  N)  ; 

case  2 

h  =  generate_targets (  Sxx(:,2)', 
N)  +  1 j *generate_targets (  Sxx(:,2)',  char (pdf_names (2 )  )  ,  N)  ; 

case  3 

h  =  generate_targets (  Sxx(:,3)', 
N)  +  1 j *generate_targets (  Sxx(:,3)',  char (pdf_names ( 3 ) )  ,  N)  ; 

case  4 

h  =  generate_targets (  Sxx(:,4)', 
N)  +  1 j *generate_targets (  Sxx(:,4)',  char (pdf_names ( 4 ) ) ,  N)  ; 

end 


char (pdf^names (1 ) ) , 
char( pdf _name  s ( 2 ) ) , 
char (pdf_names (3) ) , 
char (pdf_names (4 ) )  , 


h  =  h'  ; 


%  remove  the  mean 


realizations 


for  k5  =  1 : num_exp  %  Monte  Carlo  loop  for  number  of  noise 


number  of  loops 


loop_track  =  loop_track  +  1; 


ptheta  =  zeros  (Ntheta, 1 ) ; 

memory  for  the  prior  probabilities 

ptheta (:,1)  =  1/Ntheta; 

prior  knowledge,  we  assume  that  the 

cond_y_h  =  mp (ones (Ntheta, 1 )) ; 

variable  p(y|h) 


%  increment  on  the 

%  Vector:  Allocate 

%  RR:  initial  priors 
%  Since  we  have  no 

%  priors  are  all  equal 
%  initialize  the 


for  k6  =  1 : transmission  %  Number  of  radar  transmissions 


noise  =  sqrt ( Pn/2 ) * (randn (Ly, 1 ) + 

1 j *randn (Ly , 1 ) ) ;  %  generate  complex  noise  Pn  is  noise  power 

%  the  variance  of  the  noise  is  Pn/2 


switch  waveform_method  %  Generate  the  optimal  waveforms 
and  simulate  the  data  measurement 

case  'PWE' , 

x  input  = 

sqrt (Es ) *x_opt* sqrt (ptheta (:, k6 )) ;  %  PWE  with  proper  scaling 

y  =  conv (x_input, h)  + 

noise;  %  simulate  the  data  measurement 

case  'wideband' 


sqrt (Es ) *delta; 
noise; 


x  input  = 

%  wideband  with  proper  scaling 
y  =  conv (x_input, h)  + 

%  simulate  the  data  measurement 
end 
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%  Compute  the  "pdf" 


for  k7  =  l:Ntheta 


%  use  online  ensemble  averaging  to  compute 
%  Eg [p ( yl , . . . , yk) ] 

%  the  mp  function  is  part  of  the  multiprecision 
%  computing  toolbox 

for  k8  =  1 : num_ensemble 

%  Generate  a  target  based  on  one  of  the 
%  hypothesis 
hi 

generate_targets (  Sxx(:,k7)',  char (pdf^names ( k7 ) ) ,  N)  + 

1 j *generate_targets (  Sxx(:,k7)',  char (pdf_names (k7) ) ,  N) ; 

hi  =  h' ; 

%  compute  exp [ (y-Xg) ' KnA-l (y-Xg) ] 
yi  =  conv (x  input, hi); 

w  =  -(y  -  yi)'/Kn*(y  -  yi); 

z(k8)  =  cond_y_h ( k7 ) *exp (w) ; 

end 

%  compute  the  expected  values  of  the 
%  ensemble  average 

cond_y_h(k7)  =  mean (z) ; 

end  %  end  loop  for  computing  pdf  for  one  hypothesis 

%  the  result  is  the  pdf  for  the  estimate  of  the 
%  maximum  likelihood  detection  algorithm 
pdf  =  cond_y_h; 


%  UPDATE  THE  PRIORS 


%  RR:  update  ptheta  and  normalize  it 
ptheta ( : , k6+l )  = 

real (vpa (ptheta (:, 1) . *pdf /sum (ptheta (:, 1) . *pdf) , 2) ) ;  %  Vector  (i  x  k) 

ptheta (;, k6+l )  =  round (ptheta (:, k6+l ) *100 )/ 100 ; 
round  the  percentage  to  2  decimal  points 


end  %  end  loop  for  one  radar  transmission 


ooooooooooooooooooooooooooooooooooooooooooooooooooo 

%  the  target  will  have  the  maximum  ptheta  value 

ooooooooooooooooooooooooooooooooooooooooooooooooooo 


corresponding  to 


[pr,dec]  =  max (ptheta (:, transmission+1 ))  ; 
the  max  prior  probability 


%  dec  =  index 


ooooooooooooooooooooooooooooooooooooooooooooooooooo 


Error  calcuation 


ooooooooooooooooooooooooooooooooooooooooooooooooooo 


if  (dec  ~=  hypo) 
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err  =  err  +  1; 


end 


9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9' 

ooooooooooooooooooooooooooooooooooooooooooooooooooo 


loop  remain  = 

round ( (num_exp*target_realization*Ntheta*length (ES) *2-loop_track) ) ;  %  loop 
track 

clc; 

sprintf  ( 'loop  remain:  %d  loops.',  loop_remain)  %  display 
the  num.  of  loops  remain 

end  %  end  loop  for  one  noise  realization 

end  %  end  loop  for  one  target  realization 

error_rate ( k3 )  =  err/ (target_realization*num_exp) ;  %  compute  the 
error  rate  for  one  hypothesis 

end  %  end  loop  for  one  target  hypothesis 

error_rate_total  =  sum (error_rate) /Ntheta;  %  average  the 

error  rate  for  the  four  hypotheses 

accuracy (k2)  =  1  -  error_rate_total ;  %  compute  P(CC) 

=  1-P  (error) 

end  %  end  loop  for  one  energy  level  (Es) 

f inal_result ( kl , : )  =  accuracy;  %  assign  P(CC) 

for  the  result  of  one  energy  level 

end  %  end  loop  for  the  two  options  of  the  waveform  methods 


%  plot  the  final_result 


semilogx (ES , f inal_result, ' LineWidth' ,  2)  ; 
ylim ( [0.21]) 

xlabel ( 'Energy  per  Transmission' ) 

ylabel  (  'Probability  of  Correct  Classification') 

final_titlel  =  'Classification  Performance  with  %d  Transmissions' ; 
final_title2  =  '\n%d  Target  Realization/  %d  Noise  Realization' ; 
final_title3  =  '\n%d  Online  Target  Ensemble  Averaging' ; 

final_title  =  [ final_titlel  final_title2  f inal_title3 ] ; 

title (sprintf (final_title, transmission, target_realization, num_exp, num_ensemble) ) 
legend ( 'PWE' , ' Wideband' ) ; 


%  END  OF  CODE 


ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 
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This  code  is  written  for  generating  a  vector  of  energy  levels  specified  by  the  minimum 
and  maximum  bounds. 

function  [  Es  ]  =  energy_levels (  Es_min,  Es_max  ) 

%  create  10  energy  levels  for  1  decibel 
ES  =  zeros (1,9); 
for  i  =  1  :  9 

ES(i)  =  Es_min*i; 

end 

Es_range  =  Es_max/Es_min; 
num_decibel  =  loglO (Estrange) ; 

ES_initial  =  ES; 

%ES  current  =  zeros ( 1 , 9*num  decibel); 

ES_current  =  ES_initial; 

for  i  =  1: (num_decibel-l ) 

ES_next  =  ES_initial* 10 Ai ; 

ES_current  =  [ES_current  ES_next] ; 

end 

Es  =  ES_current; 
end 


This  is  the  code  written  to  generate  real  signals  with  both  specified  PDFs  and  specified 
PS  Vs.  This  code  was  written  by  Nichols  et  al.  [8]. 


function  [  x  ]  =  generate_targets (  Sxx,  pdf,  N) 

Sxx  (N/2  +  1 ) =0 ;  %  zero  out  the  DC  component  (remove  mean) 

Xf  =  sqrt (Sxx) ;  %  Convert  PSD  to  Fourier  amplitudes 

Xf=if ftshift (Xf ) ;  %  Put  in  Matlab  FT  format 

switch  pdf 

case  'gaussian' 

x  =  randn ( 1 , N) ; 
case  'rayleigh' 

x  =  raylrnd(10, [1,N] ) ;  %  generate  N  iid  samples  conforming  to  p (x) 
case  'exponential' 

x  =  exprnd(2, [1,N] ) ; 
case  'gamma' 

x  =  gamrnd (2,2, [ 1 , N ] ) ; 
case  'lognorm' 

x  =  lognrnd(0,  1,  [ 1 , N] ) ; 

end 

vs  =  sum ( Sxx/N) * (N/ (N-l ) ) ;  %  define  the  signal  variance  based  on  the  PSD 
x=x*sqrt (vs/var (x) ) ;  %  guarantee  new  data  match  this  variance 
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mx=mean (x) ; 


x=x-mx;  %  subtract  the  mean 

[xo, ~ ] =sort (x) ;  %  store  sorted  signal  xo  with  correct  p  (x) 


k=l ; 

indxp=zeros ( 1 , N) ;  %  initialize  counter 
while ( k) 

Rk=fft(x);  %  Compute  FT 

Rp=atan2 (imag (Rk) , real (Rk) ) ;  %  Get  phases 

x=real (ifft ( (exp (lj . *Rp) ) . *abs (Xf ) ) ) ;  %  Give  signal  correct  PSD 
[~, indx] =sort (x) ;  %  Get  rank  of  signal  with  correct  PSD 
x(indx)=xo;  %  rank  reorder  (simulate  nonlinear  transform) 
k=k+l;  %  increment  counter 
if  (indx==indxp) 
k=0  ; 

end  %  if  we  converged,  stop 

indxp=indx;  %  re-set  ordering  for  next  iter 

end 

%x=x+mx;  %  Put  back  in  the  mean 
x  =  x/sqrt (2*N) ; 

end 


This  is  the  code  written  to  generate  optimal  illumination  waveform  based  on  the  proposed 
target  hypotheses.  This  code  was  written  by  R.  Romero  [4]. 

function  [  x_opt  ]  =  opt_waveforms (  target_hypotheses  ,  time_domain) 

hh  =  target_hypotheses ;  %  reassign  the  variables 

N  =  length (hh) ; 


if  time  domain 


[~,Ntheta]  =  size (target_hypotheses ) ; 
hypotheses 

H  =  f ft (hh) ; 
target  response 

computes  the  four 

corresponding  to  the  four  target  signals 
Sf  =  abs (H) . A2 ; 

response 

Compute  the  Magnitude  squared  of  the 
targets 
else 


%  obtain  the  number  of 

%  calculate  the  fft  of  the 
%  (GAC_KNW) :  Matrix:  This 
%  target  transfer  functions 

%  (RR) : calculate  the  impulse 
%  (GAC_KNW) :  Matrix  -  We 
%  transfer  functions  of  the 
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[Ntheta,  ~]  =  size (target_hypotheses) ;  %  obtain  the  number  of 
hypotheses 

Sf  =  zeros ( length (hh) , Ntheta) ; 
for  i  =  1 : Ntheta 
H  =  hh ( i , : ) ; 

Sf(:,i)  =  if f tshift (H) ; 

end 

end 


L  =  length (Sf ) ; 

Target 

df  =  1/L; 
period) 

time  vector 
autocorrelation 

f  =  (0: (L-l ) ) *df ; 
df 

Is 

n  =  ( 0 : L-l ) ' ; 

V  =  exp  (1 j *2*pi*n*f ) ; 

kernel  trick  used  in  the  Kay  algorithm  for 
matrix  under  the  Gaussian  assumption 
by  Kay 


%  Length  of  impulse  response 
%  Length  of  the  Mag.  Sqrd  of  the 

%  Transfer  Functions 

%  delta  frequency  (frequency  sample 

%  Create  a  frequency  vector  and  a 

%  variables  needed  to  form 

%  frequency (f)  and  time (n) 

%  normalized  frequency  from  0  to  1- 

%  same  as  normalizing  sampling  time 

%  discrete  time  scale; 

%  frequency  time  vector 
%  GAC :  This  is  a  matrix  (31  x  31) 

%  This  is  a  Fourier  Transform 

%  converting  a  PSD  to  a  covariance 

%  These  are  the  eigenvectors  used 


Rh  =  zeros (L, L, Ntheta) ;  %  allocate  memory  for  the 

autocorrelation  matrix 

x_opt  =  zeros  (L, Ntheta); 


for  il  =  l:Ntheta, 
H2f  =  Sf ( : , il) ; 
one  target  hypothesis 


%  Select  one  column  out  of  the  matrix 
%  Mag.  Sqrd  of  Tranfer  function  |H(f) | A2  for 


Rh(:,:,il)  =  (V*diag (df *H2f ' ) *V' ) ;  %  form  the  autocorrelation  matrix 


[eig_max,~]  =  eigs (Rh ( : , : , il ) , 1 , ' LM' ) ;  %  find  the  maximum 

eigenvector,  which  is  the  finite  duration,  energy  constrained  x(t)  that 

%  maximizes  SNR 

x_opt(:,il)  =  eig_max/sqrt (eig_max' *eig_max) ;  %  normalize  the  energy 

in  x_opt 
end 

end 
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